AU 


rut  C'OPJ 


A  CONTINUOUS  ANALOGUE  OF  STURM  SEQUENCES 
IN  THE  CONTEXT  OF  STURM-LIOUVILLE  EQUATIONS 


L.  Greenberg 
Mathematics  Department 
University  of  Maryland 
College  Park,  Maryland  20742 


I .  Babuska 

Institute  for  Physical  Science  and  Technology 
University  of  Maryland 
College  Park,  Maryland  20742 


August  1986 
Revised  February  1988 


U  S  iiU 

SELECTE 
MAY  2  3  1988 


P 

y 

L 

DLSTRIBUnON  STATSyiNT  A  \ 

Approved  for  public  relc-aooj  ( 
Distribution  Unlimited  t 


MD88-09-LG/IB 


TR88-09 


H* 


Lnci  u; 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (WTiar i  Data  Enlarad) 


REPORT  DOCUMEHTATIOH  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

I.  REPORT  NUMBER 

MD88-09-LG/IB 

TR88-09 

1.  RCCIPICNTS  CATALOG  NUMBER 

■ 

4.  TITLE /and  Sufctlt/a) 

A  Continuous  Analogue  of  Sturm  Sequences  in 
the  Context  of  Sturm-Liouville  Equations 

J.  type  of  report  a  period  covered 

Final  life  of  the  contract 

A.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR/*) 

1  0 

L.  Greenberg-1-  and  I.  Babuska^ 

1.  CONTRACT  OR  GRANT  NUMBER/*) 

2 

ONR  N00014-85-K-01 69 

S.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Mathematics  Department  Inst-  for  Phys.  Science  6* 
University  of  Maryland  Technology 

College  Park,  MD  20742  University  of  Maryland 

College  Park,  MD  20742 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  FORK  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Department  of  the  Navy 

Office  of  Naval  Research 

Arlington.  VA  22217 

12.  REPORT  DATE 

August  1986, revised  Feb.  1988 

12.  NUMBER  OF  PAGES 

58 

14.  MONITORING  AGENCY  NAME  A  ADDRESS///  dull  rani  /real  Controlling  Olllca) 

IS.  SECURITY  CLASS,  (ol  Ihla  report) 

ISa.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (ol  Ihlm  Report) 

Approved  for  public  release:  distribution  unlimited 

17.  DISTRIBUTION  STATEMENT  (ol  (fie  obatracl  antarad  In  Bio ek  10,  II  dlllorant  from  Report) 

is.  supplementary  notes 

19-  KEY  WORDS  (Continue  on  roworoo  oldo  II  nocooaory  md  Identity  by  block  numbor) 

20.  ABSTRACT  (Continue  on  rovoroo  oldo  II  nocoooory  ond  Identity  by  block  numbor)  A  Shooting  method  IS  present  C*C 

for  finding  the  n’-U  eigenvalue  and  eigenfunction  of  a  Sturm-Lioville  equation, 
in  which  the  eigenvalue  occurs  nonlinear.lv.  The  method  is  verified  in  two 
ways:  by  applying  the  Sturm  comparision  and  oscillation  theorems  to  the 

continuous  problem;  and  by  applying  Sturm  sequences  to  a  discretization.  The 
method  works  for  general  (separated)  boundary  conditions,  and  provides  an 
a-posteriori  error  estimate  for  the  approximate  eigenvalue.  Analogues  of  the 
Sturm  comparison,  oscillat ' ■"  i  nd  separation  theorems  are  proved  for  the 
discrete  problem.  A  related  method,  which  involves  critical  lengths  in  (over) 

OD  I  j  an*71  1473  EDITION  OF  I  NOV  •»  I*  OBSOLETE 


S/N  0)02-  LF-  014-6401 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fBhan  Daia  fntarad) 


20.  the  invariant  imbedding  method,  is  shown  to  be  incorrect  for  general 
bomij.irv  conditions. 


V  V 


VI  I'lW 


A  Continuous  Analogue  of  Sturm  Sequences 
in  the  Context  of  Sturm-Liouvl 1 le  Equations 


1  x  2 

L.  Greenberg  and  I . Babuska 


Abstract . v  A  shooting  method  is  presented  for  finding  the  n  eigenvalue  and 
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parison  and  oscillation  theorems  to  the  continuous  problem;  and  by  applying 
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boundary  conditions,  and  provides  an  a-posteriorl  error  estimate  for  the 
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separation  theorems  are  proved  for  the  discrete  problem.  A  related  method. 


which  involves  critical  lengths  in  the  invariant  imbedding  method,  is  shown  to 


be  incorrect  for  general  boundary  conditions. 


Key  words:  Eigenvalue,  eigenfunction,  Sturm-Liouvi 1 le  equation,  shooting 


method,  comparison  theorems,  oscillation  theorem,  discretization,  Sturm 


sequence. 


AMS  (MOS)  classification:  65L15 


'Mathematics  Department,  University  of  Maryland,  College  Park,  MD  20742. 


"Institute  for  Physical  Science  and  Technology,  University  of  Maryland, 

College  Park,  MD  20742.  The  work  of  this  author  was  supported  in  part  by  the 
Office  of  Naval  Research,  under  contract  N00014-85-K-0169. 


£ 


«:i/  ».♦•>.* 


-f^fTT^+rf.r+rjrjrr.'w s  -jw  v.v.'v.'^v. 


1.  Sturm-Liouvil le  Equations. 


Consider  the  Sturm-Llouvi 1 le  eigenvalue  problem 


"(p(x.X)u' )'  +  q(x,A)u  =  0,  for  0  2x21, 
«0(X)u(0)  +  0  (Aju'tO)  =  0, 

ajCxjud)  +  pjixju'd)  =  o. 


(i.i) 


The  main  objective  of  this  paper  is  to  present  a  shooting  method  for  the 
eigenvalues  and  eigenfunctions  of  (1.1),  and  to  prove  its  validity.  The 
method  is  based  on  oscillation,  and  is  related  to  the  critical  (or  character¬ 
istic)  lengths  in  the  invariant  Imbedding  method.  However,  a  simple  counting 
of  critical  lengths  does  not  produce  a  correct  algorithm  (for  general,  sepa¬ 
rated  boundary  conditions).  We  shall  discuss  this  in  detail  in  §6.  Our 
method  cam  aim  for  the  n*’*1  eigenvalue  without  consideration  of  other  eigen¬ 
values.  It  provides  am  a-posteriorl  error  estimate  for  the  approximate  eigen¬ 
value.  The  method  is  a  generalization  of  that  used  by  Porter  and  Reiss  [13], 
[14],  for  the  problem 


(1.2)  •  (w-Au°(x))2 


lu(0)  =  0  =  u' (d) , 


c(x)2  (u-Au°(x))2 


u  =  0,  for  0  2  x  2  d, 


which  arises  in  acoustics.  (Here,  u  (x)  is  a  given  function. ) 

The  shooting  method  will  be  described  in  §2,  and  its  validity  will  be 
proved  using  the  Sturm  comparison  and  oscillation  theorems.  The  problem  will 
be  discretized  in  §3,  and  another  verification  of  the  shooting  method  will  be 
given  in  §4,  by  applying  Sturm  sequences  to  the  discrete  problem.  This  pro¬ 
cess  of  "talcing  the  limit"  of  a  numerical  method  (applied  to  a  discretization) 
is  sometimes  referred  to  as  finding  a  "closure  of  an  algorithm"  (see 


I 


Babuska  [1]).  Such  closures  give  Insight  into  the  numerical  method,  and  often 
lead  to  more  flexible  and  adaptable  procedures  than  the  finite  difference  or 
finite  element  method.  An  example  of  this  is  the  double  sweep  method  (see, 
for  example,  Babuska,  Prager  and  Vitasek  [4],  Babuska  and  Majer  [2],  Keller 
and  Lent  ini  [12]).  In  §5,  analogues  of  the  Sturm  comparison,  oscillation  and 
separation  theorems  will  be  given  for  the  discrete  problem.  In  §6,  we  provide 
a  cautionary  note  concerning  a  related  shooting  method  which  can  fail, 
although  its  appearance  of  validity  is  seductive. 

We  shall  need  to  make  certain  assumptions  about  the  coefficient  functions 
in  (1.1).  As  indicated  below,  these  occur  in  three  categories.  The  standard 
assumptions  will  always  be  implicitly  assumed  without  mention.  The  monotoni¬ 
city  and  limit  assumptions  occur  in  the  Sturm  comparison  and  oscillation 
theorems.  They  will  be  explicitly  assumed  when  needed.  In  the  following,  A 
will  vary  in  an  interval  (A^.Ag).  We  do  not  exclude  the  possibilities 

A1  =  ~00’  ^2  =  °°- 

Standard  Assumptions 
d 

(51)  p(x,A),  — p(x,A)  and  q(x,A)  are  continuous  functions  on 
[0,l]x(A1,A2). 

(52)  p(x,A)  >  k  >  0,  for  0  Z  x  £  1,  A1  <  A  <  A2- 

(53)  aQ(A),  PQ(A),  otj  ( A) ,  0j(A)  are  continuous  on  (Aj.Ag). 

(54)  c^U)2 + 0i(A)2  *  0,  for  i  =  0, 1  and  A1  <  A  <  A^ 

(55)  For  i  =0,1,  either  P^A)  s  0  or  P^A)  >  0,  for  A^  <  A  <  A2>  (*) 

Monotonlcltv  Assumptions 

(Ml)  For  each  x,  q(x,A)  is  a  strictly  increasing  function  of  A. 

(M2)  For  each  x,  p(x,A)  is  a  nonincreasing  function  of  A. 

(  'if  A+  =  A^  in  (L3),  then  we  assume  that  (S5)  is  valid  for  A^  S  A  <  A^ 


is  a  nondecreasing  function. 


(M3)  If  0QU)  *  0,  then  p(0,A)g|^y 

(M4)  If  P1(A)  *  0,  then  p(0,A)|±j£j  is 

We  shall  use  the  following  notation: 


nonincreasing  function. 


(1.3) 


p  (A)  =  max  p(x, A) ,  p„(A)  =  min  p(x,A), 

0£x£l  0Sx<l 

* 

q  (A)  =  max  q(x,A),  q.(A)  =  min  q(x,A). 

0<x£l  0£x<l 


Note  that  p  (A)  >  p„(A)  >  k  >  0,  by  assumption  (S2). 


Limit  Assumpt ions 


q»(A) 

(LI)  lim  -f -  =  oo. 

A-*^  p  (A) 

(L2)  lim  =  -*■ 

A-»A1  P*U) 

(L3)  There  is  a  number  A+  in  (A^.A^),  so  that  aQ(A+)0Q(A+)  £  0, 

* 

aj  (A+)/3j  (A+)  >  0,  and  q  (A+)  SO.  (If  the  coefficient  functions  in 

(  *) 

(1.1)  can  be  extended  continuously  to  A  =  A^,  we  may  take  A+  =  A^.) 

Note  that  assumption  (L3)  implies  that  for  A  =  A+,  the  operator  in  (1.1)  is 
negative  semldef inlte. 

Remark.  In  the  acoustics  problem  (1.2),  q(x,A)  is  a  decreasing  function  of 
A,  and  p(x,A)  is  an  increasing  function  of  A,  if  u  >  0,  u°(x)  >  0  and 
u>-Au°(x)  >  0.  However,  the  shooting  method  is  valid  if  the  words 
"increasing"  and  "decreasing"  are  Interchanged  in  assumptions  (Ml)  -  (M4), 
while  A^  and  A^  are  interchanged  in  assumptions  (LI)  -  (L3). 

We  shall  conclude  this  section  by  recalling  the  Sturm  comparison  and 
oscillation  theorems.  These  theorems  (with  somewhat  different  notation)  can 
be  found  in  Bocher  (5,  Chap.  3]  and  Ince  [10,  Chap.  10].  The  comparison 


If  A+  =  A^  in  (L3),  then  we  assume  that  (S5)  is  valid  for  A^  <  A  <  A 


theorems  deal  with  the  initial  value  problem  obtained  from  (1.1)  by  omitting 
the  boundary  condition  at  x  =  1.  This  can  be  formulated  as  follows: 


(1.4) 


(p(x,A),u')'  +  q(x,  A)u  =  0,  for  0  £  x  £  1, 
u(0)  =  0  (A).  u'(0)  =  -a  (A). 


Let  u(x,A)  denote  the  solution  of  (1.4).  Note  that  (for  fixed  A)  the 
zeros  of  u(x,A)  are  simple,  since  (1.4)  is  a  second  order  differential 
equation. 


Theorem  A  (First  Comparison  Theorem).  Suppose  that  (1.4)  satisfies  the 
monotonicity  assumptions  (Ml),  (M2),  and  (M3).  Then,  for  A^  <  X^,  ufx.A^) 
has  at  least  as  many  zeros  as  u(x,A^)  in  the  interval  0  <  x  £  1,  and  the 
ith  zero  of  u(x,A2)  is  less  than  the  ith  zero  of  u(x,A^). 


Theorem  B  (Second  Comparison  Theorem).  Suppose  that  (1.4)  satisfies  the 
monotonicity  assumptions  (Ml),  (M2),  and  (M3).  Let  0  <  xQ  £  1,  and  suppose 
that  u(Xq,A)  *  0  for  A^  <  A  <  A,,.  Then  p ( Xq , A ) *s  a  strictly 

decreasing  function  of  A  in  the  interval  A^  <  A  <  X^. 


Remark  1 .  The  comparison  theorems  are  usually  stated  in  terms  of  two  differ¬ 
ential  equations: 


(1.5) 


and 


(p0(x)u')'  +  q0(x)uQ  =  0, 

u0(0)  =  P0’  u6(0)  =  "V 


for  0  £  x  £  1, 


(1.6) 


(Pj(x)u')'  +  q1(x)u1  =  0,  for  0  £  x  £  1, 
UjfO)  =  0j,  u'(0)  =  -a1 , 


where  it  assumed  that 


(1.7) 


pQ(x)  £  p^x),  qQ(x)  <  (x)  and  pQ(0)j^  *  PjCO)^-. 


(For  the  second  comparison  theorem,  we  also  assume  that  Ug(x)  and  Uj(x) 
have  the  same  number  of  zeros  In  the  interval  (0,Xg),  and  Uq(Xq)  *  0, 
ul(x0)  *  0.)  However,  the  equations  (1.5)  and  (1.6)  can  be  embedded  in  a 
continuous  family  of  equations: 


(1.8) 


(p(x,A)u/)/  +  q(x,A)u  =  0,  for  0  £  x  £  1, 
u(0)  =  3(A),  u' (0)  =  -a( A ) , 


which  satisfies  the  monotonicity  assumptions  (Ml),  (M2),  and  (M3)  for 


0  <  A  <;  1,  and  such  that  p(x,0)  =  pQ(x),  q(x,0)  =  qQ(x),  p(x,l)  =  p^(x), 
q(x,D  =  qi(x),  =  jg  and  We  can  take 

p(x, A)  =  ( l-A)pQ(x)  + Ap1 (x) ,  q(x, A)  =  ( l-A)qQ(x)  + AqJ (x) ,  and  choose  a(A), 


3(A)  so  that 


Namely,  take  a(A)  =  ( l-A)pQ(O)ao01  + Ap^ (0)<*j3q.  and  3(A)  =  3Q31p(0,A).  It 
then  follows  that  the  solution  u(x,A)  of  (1.8)  satisfies 


(1.9) 


u(x,  0)  =  CqUq  (  X ) ,  u(x,  1)  =  c^jtx), 


where  cQ  =  pQ(0)31,  c1  =  p1(0)3Q.  Thus  u(x,0)  and  uQ(x)  have  the  same 

number  of  zeros  in  (0,1),  as  do  u(x,l)  and  u  (x).  Also 
u'(x0)  Ui(x) 

p(x,0) — t — ’—rr  =  p,(x) — 7—~r,  for  0  £  x  £  1.  We  shall  feel  free  to  use  either 
K  u( x, 1 )  *1  Uj(x) 

form  of  the  comparison  theorems. 


Remark  2.  If  the  strict  inequality  q^ix)  <  q^(x)  is  changed  to 

qQ(x)  i  q^(x)  in  (1.7),  then  the  conclusions  are  changed  analogously:  the 

'’irst  comparison  theorem  concludes  that  u  (x)  has  at  lea&t  as  many  zeros  as 


uQ(x)  in  the  interval  0  <  x  £  1,  and  the  i  zero  of  u^(x)  is  less  than 
or  equal  to  the  i*"*1  zero  of  Ug(x).  The  second  comparison  theorem  concludes 
u0(x0}  ulU0] 

that  pi1x,=^t  s  pi(xo,s7n^r 


Theorem  C  (Oscillation  Theorem).  Suppose  that  (1.1)  satisfies  the  monotoni¬ 
city  assumptions  (Ml)  —  ( M4 )  and  the  limit  assumption  (LI).  Then  the  eigen¬ 
values  of  (1.1)  form  an  infinite,  increasing  sequence  A  <A  ,<A  ,<..., 

m  m+l  m+ l 

which  tends  to  A^.  The  eigenfunction  y  (x),  corresponding  to  A^,  has 
exactly  k-1  zeros  in  the  interval  (0,1).  Furthermore,  suppose  that  (1.1) 
satisfies  either  the  limit  assumption  (L2)  or  (L3).  Then  the  sequence  of 
eigenvalues  begins  with  A  ,  whose  eigenfunction  <p, (x)  has  no  zeros  in 


A  Shooting  Method  for  Eigenvalues  and  Eigenfunctions 


th 


An  eigenvalue  of  (1.1)  will  be  denoted  A^,  and  called  the  k  eigen¬ 


value,  if  its  corresponding  eigenf unct Ion  has  exactly  k-1  zeros 


In  (0,1).  This  notation  has  already  been  used  In  the  Sturm  oscillation 


theorem,  which  gives  sufficient  conditions  for  A^  to  exist  for  large  k,  o 


for  all  k  >  1. 


Definition  2. 1.  For  a  given  A  in  (A^.A^),  let  u(x,A)  denote  the  solu¬ 


tion  of  the  initial  value  problem  (1.4).  Furthermore,  let 


(1) 

(2) 


N  (A)  be  the  number  of  zeros  of  u(x,A)  in  (0,1), 
o 


u(  A )  =  a  ( A  )u(  1 ,  A)  ♦f3(A)u'(l,A), 


(3) 


cr(  A )  = 


0  if  u( 1, A)u(A)  >  0  or  u(A)  =  0, 

1  if  u( 1 , A)u(A)  <  0  and  u(A)  *  0, 


(4) 


N(A)  =  N  (A)  +  <r(A) . 
o 


Remark.  In  the  above  formulation,  N ( A )  is  the  number  of  zeros  of  u(x,A) 
in  (0,1),  with  a  correction  which  depends  on  the  boundary  condition  at 


x  =  1.  If  0  (A)  *  0,  an  equivalent  formulation  is:  N(A)  =  Mq(A)  +p(A), 


where  M  (A)  is  the  number  of  zeros  of  u(x,A)  in  (0,1],  and 
o 


(A)  =  |0  if  u( 1 , A)u(A)  £  0 
ll  if  u{ 1 , A )u( A)  <  0. 


The  shooting  method  will  use  the  oscillation  counter  N(A)  in  an  essen¬ 
tial  way.  The  following  theorem  gives  the  main  properties  of  N(A). 
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Theorem  2. 1  (Shooting  Theorem).  Suppose  that  (1.1)  satisfies  the  monotoni¬ 
city  assumptions  (Ml)  -  (M4).  Let  A'  <  A"  be  numbers  in  (A^.A^)*  *. 

Then: 

(1)  The  interval  [A',  A")  contains  exactly  N(A")-N(A')  eigenvalues  of 

(1.1). 

(2)  If  N(A' )  =  j  <  k  =  N(A"),  then  the  eigenvalues  ’  ^j-t-2*  ’  ‘  '  ’ 

exist,  and  A'  <A.,<A.  „<...<  A.  <  A". 

J  +  l  j+2  k 

If  (1.1)  also  satisfies  either  the  limit  assumption  (L2)  or  (L3)  then 

(3)  For  A1  <  A  <  A^.  (1.1)  has  exactly  N(A)  eigenvalues  in  the  interval 

(ArA)(  ’. 

The  proof  of  this  theorem  is  closely  intertwined  with  the  proof  of  the 
Sturm  oscillation  theorem.  We  shall  begin  with  a  short  discussion  and  two 
lemmas,  which  will  lead  to  the  proof  of  the  theorem. 

Let  u(x,A)  be  the  solution  of  the  initial  value  problem  (1.4).  Sturm’s 
first  comparison  theorem  states  that,  as  A  increases,  u(x,A)  does  not  lose 
zeros.  It  may  acquire  new  zeros,  which  first  appear  at  the  endpoint  x  =  1, 
and  move  toward  x  =  0.  Suppose  that  the  k^  zero  appears  when  A  =  p^. 

The  p^  form  an  increasing  sequence  p^  <  pm+^  <  ...  (possibly  a  finite,  or 
even  empty  sequence).  The  solution  u(x,A)  has  exactly  k-1  zeros  in 
(0,1),  for  p^^  <  A  <  p^  (assuming  that  p^_j  and  exist).  The  second 

comparison  theorem  implies  that  p( 1 ,  A)^^  ^  is  a  strictly  decreasing  func¬ 

tion  for  p^^  <  A  <  p^-  Clearly,  this  function  decreases  from  oo  to  -oo. 

(  *) 

If  the  coefficient  functions  in  (1.1)  can  be  extended  continuously  to 
A  =  Aj ,  then  A'  may  be  taken  in  [Aj.A^h 

The  interval  [A^.A)  in  conclusion  (3)  may  be  replaced  by  (A^A)  unless 
A+  =  Aj  =  A^  in  asssumption  (L3).  See  the  remark  after  Lemma  2.2,  below. 


v>  V*  VW  V  V  V  V  V  V  V  V 

.V  AA  v -  .A/ o  ^  O  On.V/wVvrV.  vV  %*  •*  ■  '  * 


On  the  other  hand,  the  monotonicity  assumption  (M4)  states  that  p(l,X)^-^y 


oc  ( X ) 

is  a  nonincreasing  function  (assuming  that  0^(A)  *  0),  so  that  -p(l,A)^-|-^y 


is  a  nondecreasing  function.  Therefore  there  is  a  unique  number 

Vi  <  Ak <  such  that  for  A  =  p(1*A)Srfer- =  "p(1*A)lrxy> 


«j(X)u(l,X)  + PjtXJu' (1,X)  =  0.  Thus  A  =  is  the  k  eigenvalue  of 


(1.1). 


Note  that  if  £j(A)  s  0,  then  A^  =  pk  is  the  k  eigenvalue.  In 


this  case  <r(A)  is  always  zero,  and  N(A)  =  Nq(A).  Thus  N ( A )  =  k-1,  for 


“k-i  <  *  5  *v 


Lemma  2. 1.  Suppose  that  (1.1)  satisfies  the  monotonicity  assumptions 


(Ml)  -  ( M4 ) ,  and  /3  (A)  *  0.  Then  N(A )  =  k-1  for  pkl  <  A  <  A^,  and 


N(A )  =  k  for  Ar  <  A  <  pk. 


,u'  ( 1,  A) 


Proof.  Because  p(  1, A^u(i , aT~  *s  a  strictly  decreasing  function  in 


interval  pR1  <  A  <  pk>  and  -p( 1 , A)^1 1 A|  is  a  nondecreasing  function,  we 


see  that  p(l,  X)^^A^  >  -p(  1,  A)“|  | A  j-  for  pkl  <  A  <  Ak>  and 


p(1,a)u(i|aT  <  ~p(1,A)fr7XT  for  Ak  <  A  <  y  This  lmPlies  that 


|3  ( A)u(  1 ,  A)  [ax  (A)u(  1,  A)  +  (A)u' ( 1,  A)  ]  >  0  for  pkl  <  x  <  \>  and 

^(Alud.AHajfAjud.Aj  +  PjfAlu'd.A))  <0  for  Ak  <  A  <  pR.  Recalling  that 
P1(A)  >  0  and  u(A)  =  ( A)u( 1 , A)  +  ^ ( A)u' ( 1 , A ) ,  we  see  that 

u(l,A)u(A)  >  0  for  pk_^  <  A  <  Ak>  and  u(l,A)u(A)  <  0  for  Ak  <  A  c  pk_ 
Referring  to  Definition  2.1,  we  see  that  <r( A)  =  0  for  pk  <  A  <  Ak>  and 
<r( A )  =  1  for  Ak  <  A  <  pk>  Since  Nq(A)  “k-1  for  pk_^  <  A  S  ^  and 
N(A)  =  N  (A)+<r(A),  it  follows  that  N ( A )  =  k-1  for  p.  .  <  A  £  A.,  and 

O  1  K 


N(A)  =  k  for  Ak  <  A  <  pk_ 


Q.E.D. 


Lemma  2.2. 


(1)  If  (1.1)  satisfies  the  monotonicity  assumptions  (M3),  (M4)  and  the 
limit  assumption  (L2),  then  N(A)  =  0,  for  A  near  . 

(2)  If  (1.1)  satisfies  the  limit  assumption  (L3),  then  N(A+)  =  0. 


Proof.  (1)  We  shall  use  the  first  and  second  comparison  theorems  to  compare 
the  equations 


(2.1) 


(p(x.  A) ,  u' ) '  +  q(x, A)u  =  0,  for  0  <,  x  £  1, 
u(0)  =  ( A ) ,  u' (0)  =  -aQ(A). 


and 

(2.2) 

where 


(p* (A)v' ) '  +  q  ( A)v  =  0,  for  0  <.  x  <.  1, 
_y (0)  =  0(A),  v'  (0)  =  -a. 


(2-3)  a  *  p(0,A0)a0(A0),  fi{\)  =  p#(A)0o(Ao). 

Here,  AQ  is  a  fixed  number  in  (Aj.Ag),  and  A^  <  A  <  AQ.  (See  (1.3)  for 

*  * 

the  notation  p.(A)  and  q  (A).)  Note  that  p(x,A)  £  p„(A),  q(x,A)  <  q  (A), 

and  by  monotonicity  assumption  (M3), 


p(0, A) 


ao(A) 


S  P(0.Aq) 


CCq(  Ap  ) 

ihW 


P*  ( A ) 


p(0, A0)a0(A0) 
p.(A)0oUo) 


Therefore  the  comparison  theorems  apply  (as  in  Remark  2  after  Theorem  B  in 


§1).  Thus,  u(x)  does  not  have  more  zeros  than  v(x)  in  (0,1],  and 


P(1,X)UTir  *  p*U)^rrr’  u(*)  and  have  the  same  number  of  zeros 


in  (0,1)  and  u(l)  *0,  v(l)  *0.  (We  are  sometimes  suppressing  the  A, 
and  denoting  u(x,A)  =  u(x),  v(x,A)  =  v(x).) 


The  limit  assumption  (L2)  states  that  lim  9  INj-  =  Therefore,  for 

a-a2  P*UJ 


10 


3 


* 


A  near  A, ,  ■- r.T  =  -s  ,  where  s  =  s(A)  =  v  -q  , ..  v.  and  lim  s(X)  =  oo. 

1  p,  CX)  P#iX) 

Equation  (2.2)  is  equivalent  to 


(2.4) 


v"  -  s^v  =  0 

_v(0)  =  0,  v'(0)  =  -a, 


whose  solution  is 

(2.5)  v  =  8  coshsx--  sinhsx. 

s 


Recall  that  either  0Q(A)  -  0  or  /3Q(A)  >  0.  If  gg(A)  s  0,  then 

/3(A)  =  p„(A)/30(A0)  =  0,  and  a  =  p(0,A0)a(A0)  *  0.  In  this  case 

v  =  -  ^  sinhsx  has  no  zeros  in  (0,1],  so  u  doesn’t  have  any  zeros  there 

either.  If  P0(A)  >  0,  then  0  >  0.  In  fact,  /3  is  bounded  away  from  0, 

because  the  standard  assumption  (S2)  implies  that  p,(A)  £  k  >  0,  so 

/3(A)  =  p„,(A)/3q(Aq)  >  k/3Q(A0)  >  0.  Since  coshsx  >  sinhsx,  and  s— >oo  as 

A — )  A, ,  this  shows  that  v  =  0  coshsx-  -  sinhsx  >  0,  when  A  is  near 
l  s 

A,.  Again,  this  implies  that  u  has  no  zeros  in  (0,1].  Thus  we  have  shown 

i. 

that  for  A  near  A.,  N  (A)  =  0  (and  u(l,A)  *  0,  v(l,A)  *  0). 

1  o 

We  shall  now  consider  <r(A).  If  0^(A)  =  0,  then  o-(A)  =0  by  defini¬ 
tion.  Therefore  we  may  assume  (3j(A)  >  0. 


P*(A) 


v'(l,A) 
v( 1 , A) 


p*(A)s 


(0  tanh  s  -  (a/s)1 
0  -  (a/s)  tanh  sj  ‘ 


Recalling  that  p, (A)  >  k  >  0,  it  is  clear  that  p,(A)^Q^^  — »»  as 
A— »A^.  Since  p(  1 ,  A)^  |  £  p^A)^.  ,  we  also  have 

P(  1  ’ A)uTT ~~>C°  aS  t*ie  other  hand,  by  monotonicity  assumption 

(M4 ) ,  -p(l,A)|i^|  <  -p(l,A0)|l|^|,  for  A  <  AQ.  Therefore,  for  A  near 

V  P,I’>liTT7#  >  -p(1-A)lrTXT-  8nd 

u( 1 , A)u(A)  =  u(l,A)[a  (A)u(l,A)  +  0j (A)u' ( 1,  A) ]  >  0. 
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This  shows  that  <r(A)  =  0  for  A  near  A^.  Therefore 

N(A)  =  N  (A)+<r(A)  =  0,  for  X  near  A,. 

o  1 

(2)  We  shall  compare  the  equations 


(2.6) 


and 


(p(x, A+),u' )'  +  q(x,A+)u  =  0,  for  0  £  x  £  1, 
u(0)  =  30(A+),  u'(0)  =  -aQ(X+). 


(2.7) 


*(P*(X+)v'  )'  +  q  (X+)v  =  0,  for  0  5x51, 
v(0)  =  3(A),  v' (0)  =  -a, 


where 


(2.8) 


a  *  p(0, X+)aQ(X+),  3  =  p,(X+)P0(X+). 

Since  p(x,X+)  £  p»(A+),  q(x,X+)  5  q  (A+)  and  p(0,X+)^~|  =  p»(X+)|,  the 

* 

first  and  second  comparison  theorems  apply.  By  assumption  (L3),  —  5  0. 


Therefore  q  /**!  =  -s2,  where  s  *  V^-q> 


P* 


=  -s  ,  where  s  -  v  510(1  equation  (2.7)  is  equlva 


lent  to 


(2.9) 


v"  -  s2v  =  0 

v(0)  =  3,  v' (0)  =  -a, 


whose  solution  is 


(2.10)  v  =  3  cosh  sx  -  ^  sinhsx,  if  s  >  0, 

or 


(2. 11)  v  =  -ax  +  3.  if  s  =  0. 

By  assumption  (L3),  ag(A+)3g(A+)  5  0,  which  implies  that  a3  5  0.  (Of 


course,  a  and  3  cannot  both  be  zero.)  Therefore  (2.10)  and  (2.11)  show 
that  v(x)  has  no  zeros  in  (0,1).  Consequently  u(x)  has  no  zeros  in 


(0,1],  and  N  (X  )  =  0. 

o  + 


If  s  >  0,  then  pBVJtTT  =  p.(A 


♦KB 


tanh  s  -  (a/s) 
-  (a/s)  tanh  s 


]  >  o, 


since 


ofi  *  0.  If  s  =  0,  then  =  P.U+)[r^]  *  0. 


Since 


u'  ( 1 1  V' m  u'  ( 1 ) 

p(  l,X+)^n~J~  -  P,(^)^n  >  this  shows  that  p(  1,  A^I^yy-  £  0.  On  the  other 


hand,  a^  (A+)/3j  (A+)  >0  by  assumption  (L3).  This  implies 

-'’‘••VgrricT 5  °'  30  that  2  -ptl- *nd 

u(l,X+)u(X+)  =  u(l,X+)[«1(A+]u(l,X+)  +  Pj(X+)u'  (1,A+)]  £  0.  Noting  that 


u(l,A+)  *  0,  this  shows  that  <r(A+)  =  0.  Therefore 


N(X+)  =  N0(X+)  +  <r(X+)  =  0. 


Q.E.D. 


Remark.  By  the  previous  lemma,  if  (1.1)  satisfies  the  limit  assumption  (L3), 


then  N(A+)  =  0.  If  (1.1)  also  satisfies  the  monotonicity  assumptions 


(Ml)  -  (M4) ,  and  A1  <  A+,  then  N(A)  =0  for  <  A  <  A+.  This  will 


follow  from  Theorem  2.1  (1),  which  implies  that  N(A)  is  a  nondecreasing 


function.  Also,  if  A^  <  A^  then  N(A)  =  0  for  A^  <  A  <  A^ .  This  will 


follow  from  Theorem  2.1  (3).  However,  it  can  happen  that  A+  =  A^  =  A^ . 


this  case,  N(A)  =  1  for  A  near  A^.  An  example  of  this  is  the  Sturm- 


Liouville  problem 


(2. 12) 


u"  +  A  u  =  0,  for  0  <  x  <  1, 


u' (0)  =  0  =  u' (1), 


where  A^=0<A<oo=  A In  this  case,  the  first  eigenvalue  is 


Aj  =  0  =  Aj,  with  eigenfunction  ^(x)  =  1. 


Proof  of  Theorem  2.1. 


(1)  Suppose  that  f3j(A)  *  0.  THen  Lemma  2.1  implies  that  N(A)  is  a 


piecewise  constant  function  with  jump  discontinuities  at  the  points 


A,  :  N(A,_+e)  =  N ( A,  )  +  1,  N(A,_-e)  =  N(A,J.  Therefore,  for  A'  <  A"  in 


N 


i 


(A^.A^),  N(A")-N(A')  equals  the  number  of  jump  discontinuities  of  N(A)  in 
[A', A"),  which  equals  the  number  of  eigenvalues  in  [A', A").  If  0^(A)  a  0, 
then  the  same  is  true,  but  the  eigenvalues  are  A^  = 

(2)  Suppose  N(A'  )  =  J  <  k  =  N(A").  By  (1),  there  exist  exactly  k-j 
eigenvalues  in  the  interval  [A', A").  If  0^(A)  a  0,  then  Lemma  2.1  Implies 
that  A'  <  Aj+1  <  A^  <  A".  The  same  is  true  if  0^(A)  a  0,  where  the  eigen¬ 
values  are  A^  =  Thus  A'  £  A,+^  <  *j+2  <  <  <  •  ( If  J +  1  =  k, 


this  sequence  contains  only  one  eigenvalue. ) 

(3)  Lemma  2.2  implies  that  there  exists  AQ,  A^  <  A^  <  A g,  such  that 
N(Aq)  =  0.  (We  may  assume  A^  <  A^,  unless  A+  =  A^  =  A^.  In  this  latter 
case,  conclusion  (1)  remains  valid  if  A'  =  A^. )  If  A^  <  Aq,  then  (1) 
implies  that  there  are  no  eigenvalues  in  the  Interval  (A^.A^),  and  N(A)  =  0 
for  A^  <  A  <  Aq.  Thus  (3)  is  true  for  Aj  <  A  S  A^. 

We  now  drop  the  assumption  A1  <  XQ  and  assume  AQ  <  A.  Any  eigenvalues 
in  [ Aj , A)  are  contained  in  [AQ,A),  and  the  number  of  these  eigenvalues  is 


N ( A )  -  N(Aq)  =  N(A). 


Q.E.D. 


Remark.  The  Sturm  oscillation  theorem  implies  that  if  (1.1)  satisfies  the 
monotonicity  assumptions  (Ml)  -  (M4)  and  the  limit  assumption  (LI),  then 
lim  N(A)  =  oo. 


We  can  now  describe  the  shooting  method. 

STEP  0.  Find  values  L^  <  R^,  such  that  N(Lq)  =  n-1  and  N(Rq)  =  n. 

This  implies  that  L„  £  A  <  R_,. 

K  0  n  0 

STEP  k.  For  given  values  L^^  <  R^_j,  with  N(L^_j)  =  n-1,  N(R^_j)  =  n, 
find  values  L^.R^.  such  that  N(L^)  =  n-1,  N(R^)  =  n. 

h.-i  s  Lk  <  "k  5  "k-r  *nd  "k '  Lk  <  \-i  ‘  Lk-r 


STOP  when  <  x,  where  x  is  a  given  tolerance. 


To  implement  the  above  steps,  we  need  an  initial  value  solver  to  compute 

the  solution  u(x,A)  of  (1.4),  when  A  =  L^.R^.  This  allows  us  to  calculate 

N(L^)  and  N(R^).  We  also  need  a  nonlinear  solver  to  calculate  the  new 

values  Ljt+j,Rjc+^.  We  have  in  mind  a  combination  of  the  bisection  method  for 

the  integer  function  N(A),  and  some  other  method  (such  as  the  secant  method) 

for  the  continuous  function  u(A)  =  (A)u( 1 , A)  + 0^ ( A)u' ( 1 , A) .  Denoting 

=  (L^+R^)/2,  the  bisection  method  would  set  =  L^,  R^+j  =  if 

NlMj^)  =  n,  or  L^+1  =  M^,  R^+1  =  R^  if  N(M^)  =  n  -  1.  The  other  part  of 

the  solver  would  be  an  iterative  method  to  solve  u(A)  =  0.  The  approximate 

eigenvalue  A^  will  be  either  the  midpoint  of  the  last  interval  [L^.R^],  or 

the  last  approximation  A^  found  by  the  iterative  method.  The  approximate 

eigenfunction  #>n(x)  will  be  the  solution  of  (1.4),  with  A  =  A^.  Usually, 

this  will  have  been  calculated  already,  and  no  further  work  will  be  needed. 

STEP  0  can  be  carried  out  either  by  using  estimates  of  A  which  the  user 

n 

might  have,  or  by  using  a  related  boundary  value  problem  whose  eigenvalues  are 
known. 

Assuming  that  the  initial  value  problems  are  solved  exactly,  the  method 
gives  a  sharp  a-posteriori  error  estimate  for  the  eigenvalue  A^.  Of  course, 
the  differential  equations  will  be  solved  numerically.  An  effective  implemen¬ 
tation  of  this  method  must  relate  the  accuracy  of  the  initial  value  solver 
(governed  by  an  input  tolerance  parameter)  to  the  value  R^  -  L^,  and  to  the 
accuracy  of  the  nonlinear  solver  for  finding  L^,R^. 

If  we  are  only  interested  in  the  eigenvalue  A^,  and  not  the  eigenfunc¬ 
tion  ^(x),  then  we  need  only  concern  ourselves  with  the  count  of  the  zeros 
of  u(x,A)  and  the  correction  term  <r.  This  can  be  obtained  by  solving 


,  JV  i  v~_  w  -r„  •r,  ■*%  •>' «,'•  *  ,  *  .  *•  *  •  _  -  _  n  «  n  •  "i  *  « 


various  transformed  formulations  (such  as  that  used  in  the  invariant  imbedding 
method) . 


The  method  resembles  a  count  of  the  critical  lengths  in  the  invariant 
imbedding  method  (see  Scott  (15,  Chap.  5]).  However,  a  simple  count  of  the 
critical  lengths  (with  no  correction  term)  does  not  produce  a  correct  algo¬ 
rithm,  for  general  boundary  conditions.  We  shall  return  to  this  point  in  §6 


Remark.  We  can  obtain  another  version  of  the  shooting  theorem  if  the  words 
"increasing"  and  "decreasing"  are  interchanged  in  the  monotonicity  assumptions 
(Ml)  -  (M4) ,  while  A1  and  A^  are  interchanged  in  the  limit  assumptions 
(LI)  -  (L3).  In  this  case,  the  Sturm  oscillation  theorem  would  declare  the 
existence  of  a  decreasing  sequence  of  eigenvalues  A^  >  >  ^m+2  >  • • • » 


which  tends  to  A^.  There  is  no  change  in  the  definition  of  N(A).  The  con¬ 


clusions  in  the  shooting  theorem  would  be  changed  to  the  following: 

(1)  The  interval  (A', A"]  contains  exactly  N(A')-N(A")  eigenvalues 
of  (1.1). 


(2)  If  N(A' )  =  j  >  k  =  N(A") ,  then  the  eigenvalues  *k+2 . 


exist,  and  A'  <  A^  <  A^_^  <  ...  <  A^+2  <  5  *"• 


(3)  For  Aj  <  A  <  A 2>  (1.1)  has  exactly  N(A)  eigenvalues  in  the 


interval  (A.AgJ. 

This  version  of  the  shooting  theorem  would  apply  to  a  Sturm-Liouvi 1 le 
equation  such  as  the  acoustics  problem  (1.2).  This  new  version  of  the  theorem 
easily  follows  from  the  old  version  by  considering  the  functions 


p(x,A)  =  p(x,-A),  q(x,  A)  =  q(x,-A),  a^A)  =  a^-A)  and  3j(A)  =  p^(-A). 
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3.  The  Discrete  Problem. 


In  this  section,  we  shall  discretize  the  boundary  value  problem  (1.1). 

We  shall  generate  a  difference  scheme  by  using  finite  elements.  This  is  a 

2 

convenient  method  which  guarantees  a  local  0(h  )  error  and  a  symmetric  dif¬ 
ference  matrix.  In  §4,  we  shall  apply  Sturm  sequences  to  the  discrete 
problem.  This  will  lead  to  another  proof  of  the  shooting  theorem,  and  will 
account  for  the  similarity  between  the  shooting  method  and  the  Sturm  sequence 
algorithm. 

Recall  that  the  energy  inner  product  B(u,v)  =  B(A;u,v)  for  the  problem 

(1.1)  is  given  by 

(3.1)  B(u,v)  =  -  gj|£j.p(0,A)u(0)v(0)  +  gj£j-p(l,A)u(l)v(l) 


,1 


(pu' v'-quv)dx 


(•) 


A  weak  solution  of  (1.1)  (for  a  fixed  A)  is  a  function  u  in  the  Sobolev 
space  H*[0,1],  such  that 

(3.2)  B(u,v)  =  0,  for  all  vcH^O.l]. 

If  (3.2)  admits  a  nontrivial  solution  Uq(x)  for  a  particular  value  A  =  A^, 

then  Aq  is  an  eigenvalue,  and  Uq(x)  a  corresponding  eigenfunction.  (If 

B_( A)  s  0,  then  ^  is  set  equal  to  0  In  (3.1),  and  H*[0,1]  is  replaced 
0  Po 

by  the  subspace  of  functions  v  e  H*[0,1),  such  that  v(0)  =  0.  The  case 
0j(A)  =  0  is  treated  similarly.  We  shall  carry  out  the  calculations  in  the 
generic  case  BQ(A)  *  0>  0^(A)  *  0.) 

The  problem  will  be  discretized  using  piecewise  linear  functions,  with 

(  'strictly  speaking,  the  energy  inner  product  for  the  operator  in  (1.1)  is 
the  negative  of  the  inner  product  (3.1).  But  this  is  irrelevant  for  the 
equation  B(u,v)  =  0. 


uniform  mesh  h  =  1/n.  Consider  the  partition  x_  <  x,  <  . . .  <  x  of  the 

interval  1 0 , 1 ]  by  the  nodes  x^  =  ih.  The  finite  element  spac6  Sh  is  the 

space  of  continuous  functions  on  [0,1]  which  are  linear  on  each  Interval 

[Xj_^,Xj].  The  inner  product  (3.1)  will  now  be  restricted  to  S^,  and  the 

integrals  in  (3.1)  will  be  approximated  by  quadrature  formulas.  We  shall  use 

r1 

the  midpoint  rule  for  the  integral  pu'v'dx,  and  the  trapezoid  rule  for  the 

r1 

integral  quvdx.  This  defines  an  inner  product  B,  (u,v)  on  S,  .  The 
J0  h  h 

finite  element  solution  is  a  function  u  €  S,  ,  such  that 

h 


(3.3) 


B.  (u,v)  =  0  for  all  v  e  S.  . 
h  h 


A  basis  v,,,v . ,v  for  S.  can  be  obtained  as  follows: 

0  1  n  h 


(3.4a) 


v0(x)  = 


f-x/h+  1,  for  Xq  S  x  £  Xj, 
D  el sewhere . 


(3.4b) 


v^x)  = 


(x-x^J/h  +  1,  for  x^_j  5  x  i  Xj, 
-(x-Xj)/h  +  1,  for  Xj  S  x  S  xi+l’ 
0  elsewhere, 


for  1  S  i  S  n  -  1 . 


(3.4c) 


v  (x)  = 
n 


(x-x  )/h+l,  for  x  ,  £  x  £  x  , 
n  n-1  n 

0  elsewhere. 


This  basis  is  uniquely  determined  by  the  property 


(3. 4d) 

It  also  satisfies 


v^(Xj)  =  5^ j  (0  <  i.J  <  n). 


(3. 4e) 


The  support  of  v^(x)  is  contained  in  the  one  or  two 


intervals  which  contain  x 


i' 


n 


A  function  u  6  Sh  can  be  expressed  in  the  form  u(x)  =  £  u^v^(x),  where 


i=0 
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Uj  =  u(Xj).  The  equation  (3.3)  is  now  equivalent  to  the  system  of  equations 


n 


(3.5) 


^uiBh(v1,Vj)  =0,  J  =  0,1 . n. 


1=0 


Because  of  property  (3.4e),  ®h^Vi’VJ^  =  ®  I  i  —  J I  >  1  -  A  simple  calcu¬ 

lation  shows  that 


(3.6) 


Pi 

r ,  ,  s  OLq  ^  2  h 

Bh(v0,v0)  "  00  P0  h  ~  2  q0’ 

P  _i 

Bh(vn’  vn>  '  r,  Pn  *  “h1  '  1  V 


Bh(vl,vl’  =  K  tp,  1  +  P  1515  n-1. 

’  2  ‘*5 

p  1 

ywi'-T' 


Here,  p^  =  Pj(A)  denotes  p(Jh,A),  where  J  is  an  integer  or  half-integer, 
and  similarly  for  q^  =  q^(A).  Multiplying  the  equations  (3.5)  by  h,  we 
obtain  a  system  of  equations: 


(3.7) 


AU)u  =  0, 


where  u  =  (u„,u . ,u  )  and  A(A)  is  the  matrix 

0  1  n 


(3.8)  AU)  = 


(bo~aoJ 


-b„ 


■bo  (Vbrai) 


(bl+b2  a2)  "b2 


-b  „  (b  _+b  ,  -a  ,)  -b  . 

n-2  n-2  n-1  n-1  n-1 


n-1 


(b  .-a  ) 
n-1  n 


The  matrix  coefficients  here  are 


(3.9) 


aCr 


!r<o(A)  +Wpou)’ 


a  (X) 
n 


r  V*1  -  Irm  »nU)’ 


ai  (X) 


h  q^(X),  for  1  S  i  £  n-1. 


b1  (X) 


p  i(X),  for  0  S  i  S  n-1. 
l  ♦- 
2 


Recall  that  we  have  assumed  (for  i  =0,1)  that  either  ^(X)  =  0,  or 


j8^  t  X )  >  0.  If  ( X )  a  0,  then  u^  =  0  and  the  first  row  and  column  In 


A ( X )  are  omitted.  If  /3,(X)  a  0,  then  u  =0  and  the  last  row  and  column 

1  n 


are  omitted.  Thus  A(X)  is  an  mxm  matrix,  where  m  can  be  n-  1,  n  or 
n+1.  We  shall  continue  to  confine  our  calculations  to  the  generic  case 


/3^(X)  >  0,  for  1  =0,1.  In  this  case,  A(X)  is  an  (n+l)x(n+l)  matrix. 


The  standard  assumptions  (SI)  -  (S3)  in  §1  imply  that  the  a^X)  and 


bj(X)  are  continuous  functions,  and  bj(X)  ^  k  >  0.  The  monotonicity 


assumptions  (Ml)  -  (M4)  imply  that  the  functions  a^(X)  are  strictly  increas¬ 


ing  and  the  bj(X)  are  non increasing. 


Remark.  The  finite  difference  method  leads  to  almost  the  same  difference 


scheme  as  (3.7).  Using  the  difference  operator  Au,  =  (u  j-u  j)/h, 

1  i+-  i-- 

2  2 


and  approximating  the  differential  equation  in  (1.1)  by  the  difference  equa¬ 
tions  A(p^Au. )  + qjuj  =  0.  we  obtain  all  of  the  equations  in  (3.7)  except 
the  first  and  last  equations  (corresponding  to  the  first  and  last  rows  of 
A(X)).  For  the  discretization  of  the  boundary  conditions  in  (1.1),  a  standard 


procedure  is  to  discretize  the  derivative  u'(0)  by  (u^-u  ^)/2h,  the  deri¬ 
vative  u'(l)  by  (u  -u  j)/2h,  and  to  extend  p(x,X)  slightly  outside  of 
the  interval  [0,1]  by  reflecting  values  about  the  endpoints  x  =  0  and 

x  =  1.  If  we  carry  out  this  procedure,  we  obtain  new  first  and  last  equations 

2 

which  differ  by  0(h  )  terms  from  the  old  first  and  last  equations  in  (3.7). 


'  •  '  t.  »-« >rw  >r,  -j>  .■  v  w>  '> V'-'V- .v 


Therefore,  we  may  regard  the  first  and  last  equations  in  (3.7)  as  discretiza¬ 
tions  of  the  boundary  conditions  in  (1.1).  We  may  also  regard  the  first  equa¬ 
tion  in  (3.7)  together  with  the  equation  uQ  =  0q(A),  as  a  discretization  of 
the  initial  conditions  in  (1.4). 

If  we  set  v  =  pu' ,  then  the  second  order  equation  (1.4)  is  converted  to 
a  first  order  system 


(3. 10) 


u'  =  v/p, 

•  v'  =  -qu, 

u(0)  =  $0,  v(0)  =  -p(0)aQ. 


The  first  n  equations  (3.7)  are  equivalent  to  the  difference  scheme 


(3. 11) 


u,  =  u 


v. . 


1  “1-1  p  i  1 


i  -; 


vi  ■  Vi  '  hqi-iVr 


where 


-  V 


This  is  an  implicit,  general  one-step  method, 


which  converges  to  the  solution  of  (3.10).  (See,  for  example,  Halrer,  Norsett 


and  Wanner  [9].)  Therefore  the  discrete  solution  u^  of  (3.7)  (with  the  last 


equation  omitted)  converges  to  the  continuous  solution  u(x,\)  of  (1.4),  and 


the  discrete  derivative  u'  =  — —  Ul~1 


converges  to  u'(x,h). 


4.  Application  of  Sturm  Sequences. 


Let  A ( A )  be  the  finite  element  matrix  (3.8).  The  Sturm  sequence 


SQ (A),S^(A) . Sn+^(^)  for  A(A)  is  defined  as  follows.  Sq(A)  =  1,  and 


for  1  <  i  <  n+1,  S^(A)  is  the  i 


th 


principal  minor  of  A(A): 


(4. 1)  Si(A)  = 


(b0  ao)  -b0 

-bo  (VVV  ~bi 


-bl  (b1+b2-a2)  -b2 


“bi-3  (bi-3+bi-2*al-2)  ~bl-2 


'1-2 


(bi-2+bi-rai-r 


(In  the  case  i  =  n+ 1,  b.  =  b  is  defined  to  be  zero. )  The  Sturm 

1-1  n 

sequence  should  not  be  computed  directly  from  its  definition  (4.1),  but  from 
the  recursion  relation 


(4.2) 


si+i  =  (bi-i+Vai)si  "  bi-isi-r  for  1  *  1  *  n~1- 


S  .  =  (b  -a  )S  -  b  ,S 
n+1  n-1  n  n  n-1  n-1 


which  follows  from  the  cofactor  expansion  of  the  determinant  by  its 

last  row.  Note  that  S  ,  =  det(A(A)l.  The  main  significance  of  the  Sturm 

n+1 

sequence  is  that  it  can  be  used  to  approximate  the  zeros  of  de t [ A ( A ) ) 

(which,  in  turn,  approximate  the  eigenvalues  of  (1.1)).  See  Stoer,  Bullrsch 
(17)  for  the  application  of  Sturm  sequences  to  linear  eigenvalue  problems. 


Definition  4.1.  For  a  given  A  in  (A^.Ag),  c(A)  is  the  number  of  sign 

changes  in  the  sequence  Sq(A),Sj(A) . S^^fA),  after  the  zero  terms  (if 

any)  have  been  omitted  from  the  sequence. 


,li 


The  following  theorem,  which  Is  proved  in  Greenberg  [7,  §4],  gives  the  main 
facts  about  the  Sturm  sequence  (4.1). 


Theorem  4.1.  Let  A(A)  be  the  finite  element  matrix  in  (3.8).  Suppose  that 
(1.1)  satisfies  the  monotonicity  assumptions  (Ml)  -  (M4).  Then: 

(1)  For  any  numbers  A'  <  A"  in  (A^.A^),  det[A(A)]  has  exactly 

c(A")-c(A')  different  zeros  in  the  interval  [A', A"). 

(2)  The  zeros  of  S^(A)  and  S^_^(A)  are  interlaced,  for  1  <  i  <  n+1. 

If  (1.1)  also  satisfies  either  the  limit  assumption  (L2)  or  (L3),  then: 

(3)  For  any  AQ  in  (A^A,,),  det[A(A)]  has  exactly  c(AQ)  different 

(  *) 

zeros  in  the  interval  (A^.A^) 

If  (1.1)  satisfies  assumptions  (Ml)  -  (M4),  (LI)  and  either  (L2)  or  (L3), 


then: 


(4)  det  [  A(  A)  ]  has  exactly  n+1  different  zeros  in  (Aj.A,,) 


We  shall  now  use  Sturm  sequences  to  study  the  discrete  problem.  The 
initial  value  problem  (1.4)  is  discretized  by  the  first  n  equations  of  the 
system  (3.7),  together  with  the  initial  condition  u^  =  f?o(A).  Let  u^ 

(0  <  i  <  n)  denote  the  solution  of  the  discrete  initial  value  problem. 


Theorem  4.2.  u.  =  t^S^  (for  0  <,  i  <  n),  where  t^  = 


(  *If  A+  =  Aj  in  limit  assumption  (L3),  then  the  interval  (A^.Aq)  should 
be  replaced  by  (Aj.Aq)  in  conclusion  (3),  and  ( A} , A^ )  replaced  by 
[A  .Ag)  in  conclusion  (4). 

1  ^his  assumes  that  (3q(A)  >  0.  See  the  remark  following  the  theorem,  for 
the  case  pg ( A )  =  0. 


1 


V 


Proof.  Consider  the  first  (i+1)  equations  of  (3.7): 


■(Wuo 


■bouo  *  (VWU1 


(4.3) 


'bi-2Ui-2  +  tbl-2+bl-l~al-l)ul-l  “  bl-lUi 


'bi-iui-i  *  (bi-i*brai)ui 


Vi.r 


We  can  solve  for  u^  using  Cramer’s  rule,  to  obtain  u^ 


biui»isi 


,  which 


can  be  written: 


(4.4) 


Vl  1  Ui 


SU1  bi  Si' 


First  suppose  that  all  Sj  *  0  and  all  u^  *  0.  Letting  u^/Sj  =  tj. 


(4.4)  shows  that  tJ  +  J  =  t  ^/b^ .  Thus  tJ  =  tjj/bj^  =  ti-2/^bi-lbi-2^  ~  ••• 


=  t0/^bi-lbi-2‘ - ‘ b0^ ‘  Slnce  uq  =  ^0  and  S0  =  t0  ^0  and 


h  *  V(bl-lbl-2  -  b0>' 


Now  consider  the  case  where  some  u^  or  Sj  are  zero.  Note  that  two 


consecutive  u^  cannot  be  zero,  and  two  consecutive  Sj  cannot  be  zero. 


This  follows  from  the  recursion  relations  (4.2)  and 


(4.5) 


-bi-iVi*<bi-i*brVui-biVi =  °- 


u1  =  0  =  u1+j,  then  (4.5)  Implies  that  u^_j  =  (Recall  that 


b,  j  =  p  i  *  0).  Similarly,  0  =  uJ_1  =  u^2 
'*£ 


.  =  Uq,  which  contradicts 


uQ  =  *  0.  In  the  same  way,  (4.2)  implies  that  two  consecutive  are  not 


zero.  Equation  (4.4)  may  be  written  ujSj  +  i  =  b^u^S^,  which  shows  that 
u^  =  0  if  and  only  if  =  0.  If  u^  =  Sj  =  0,  then  equations  (4.2)  and 


(4.5)  imply 


CW8 


(4.6) 


1-1..  q  _  h2  Q 

bt  i-r  si+i  “  bi-isi-r 


Therefore  =  r— ^ so  that  t.  ,  =  t  /(b.b  ),  even  if 

Si+i  ^1^1-1  Sj-i  1  +  1  1-1  1  1-1 


ui  =  si  -  °- 


Q.E.D. 


Remark  1.  In  Theorem  4.2,  we  have  assumed  that  £q(A)  >  0-  If  0Q(A)  3  0. 
then  the  first  equation  in  (3.7)  is  omitted,  uQ  =  0,  and 


<VbrV  -bi 

-bj  (b^b2-a2)  -b2 


(4.7)  S  - 


-bi-2  (bl-2*bl-rai-l)  -bi-i 

-bi-i  lbi-i*brai) 


In  this  case  (4.4)  is  replaced  by 


(4.8) 


Ui+i  _  1_  _ui_ 
Si  '  bt  Si.j’ 


and  so  Uj,  =  t^S^^  (for  1  <,  i  S  n),  where  t^  =  t^/lb^h^  ^.  .  .  b} ) ,  and 

t  =  Uj  *  0.  Thus,  if  0Q  =  0,  then  all  t1  have  the  same  sign  as  u^.  If 

/3g  >  0,  then  Theorem  4.2  shows  that  all  t^  >  0.  In  both  cases,  the  solution 

(Uq.u^ . u^)  of  the  discrete  initial  value  problem  has  the  same  number  of 

sign  changes  as  the  Sturm  sequence  S^.S^ . S^.  (Note  that  if  0Q  >  0, 

there  is  an  additional  function  S  ,,  which  has  not  been  included  in  the 

n+I 

last  sequence. ) 

Remark  2.  Theorem  4.2  indicates  that  we  can  expect  overflow  to  occur,  if  we 
attempt  to  calculate  for  small  h.  For  u^  can  be  expected  to  have 

moderate  size,  since  it  is  the  discrete  solution  of  the  initial  value  problem 


raws 


,«4_ 


’l  ~  W 


fl  1 

r  i  i 

b0bl  -bl-l  = 

[poj 

£np(x)dx 

fl  1  hJ0 

Ur 


»  |^-Je  •  This  observation  illustrates  how  how  the  closure  method 

can  give  insight  into  a  numerical  process. 

Note  that  the  matrix  A(A)  and  all  the  functions  derived  from  it  depend 
on  the  mesh  size  h.  We  shall  now  exhibit  this  dependence  explicitly: 


AC  A)  =  A(h,  A) ,  S1  =  Sj  (h,  A) ,  c(A)  =  c(h,A),  Ul  =  u^h.A). 


Theorem  4. 3.  Suppose  that  (1.1)  satisfies  the  monotonicity  assumptions 

(Ml)  -  (M4),  Aj  <  A  <  A^  *,  and  A  is  not  an  eigenvalue  of  (1.1).  Then 

lim  c(h, A)  =  N(A) . 
h-+0 


Proof.  We  shall  assume  that  /3Q(A)  >  0-  The  case  /3q(A)  s  0  is  similar. 

(See  the  remark  at  the  end  of  the  proof. ) 

Let  co(h,A)  be  the  number  of  sign  changes  in  the  sequence 

S  (h, A) , S. (h, A) , . S  (h, A) .  (Note  that  the  last  function  S  .  has  been 

U  1  n  n+1 

omitted.)  By  Theorem  4.2,  CQ(h,A)  equals  the  number  of  sign  changes  in  the 

sequence  u_,u . u  (because  all  t.  >  0).  If  u(l,A)  *  0,  then,  as 

u  l  n  l 

h— >0,  CQ(h,A)  approaches  the  number  of  zeros  of  u(x,A)  in  (0,1).  Thus 

11m  c  (h,A)  =  N  (A).  (The  possibility  u(l,A)  =  0  will  be  considered  below.) 

.  °  o 

h-»0 

We  shall  now  look  for  the  relation  between  the  last  function  S  ,  and 

n+1 

the  correction  term  <r(A).  By  the  recursion  relation  (4.2), 

2 

S  ,  =  (b  -a  )S  -b  ,S  Using  the  relations  u.  =  t.S,  and 

n+1  n-1  n  n  n-1  n-1  6  1  i  i 

tj  +  j  *  t^/b^  from  Theorem  4.2,  we  obtain 


If  the  coefficient  functions  in  (1.1)  can  be  extended  continuously  to 

A  =  A, ,  then  Theorem  4.3  is  valid  for  A,  <,  A  <  A„. 

1  12 


yr  v*  ■„%  --<  vvl'.'w.vi-o.''o  v  * '  - ,- 


Recall  that  u(A)  =  otj(A)u(l,A)  +  0j(A)u' (1,  A).  Clearly  a^un  +  pjU^— ■ »u 
as  h — >0.  Since  we  have  assumed  that  A  Is  not  an  eigenvalue  of  (1.1), 
u(A)  *  0.  Equation  (4.9)  shows  that  sgn  =  sgn  u(A),  for  small  h.  We 

shall  now  consider  two  cases,  according  to  the  possibilities  that  u(l,A)  =  0 
or  not. 

Case  1 .  u(l,A)  *  0. 

For  small  h,  cQ(h,A)  =  Nq(A),  sgn  Sn+^  =  sgn  u(A)  and  sgn  =  sgn  * 

sgn  u(l,A).  Thus  sgn  SnS^+1  =  sgn  u(l,A)u(A).  Referring  to  Definition  2.1, 

this  shows  that  <r(A)  =0  if  S  S  ,  >  0,  and  <r(A)  =  1  if  S  S  <0. 

n  n+1  n  n+1 

Therefore  c(h,A)  =  N(A),  for  small  h. 


27 


In  this  case,  <r(A)  =  1  and  N(A)  =  N  (A)  +  1.  Since  u(l,A)  =  0, 

o 

u'(l,A)  *  0.  We  shall  suppose  that  u'(l,A)  >  0.  (The  case  u'(l,A)  <  0  is 
similar. )  Since  u(A)  =  (A)u' ( 1, A) ,  and  u(A)  *  0,  it  follows  that 

^(A)  >  0  and  u(A)  >  0. 

For  small  h,  the  approximate  solution  (u_,u . u  )  will  have  N  (A) 

0  1  n  o 

sign  changes  corresponding  to  the  interior  zeros  of  u(x,A).  (Since 
u(l,A)  =  0,  there  may  also  be  a  spurious  sign  change  near  the  end  of  the 
sequence. 1  The  Nq(A)  sign  changes  will  occur  in  an  initial  subsequence 
Uq, Uj, . . . , u^ ,  where  the  node  ih  is  near  the  largest  interior  zero  of 
u(x, A) . 

Since  u'(l,A)  >  0  and  u(l,A)  =  0,  there  is  a  point  xQ  in  (0,1), 
such  that  u(x,A)  <  0  for  xQ  S  x  <  1,  and  u' (x,A)  >  0  for  xQ  <  x  <  1. 
Correspondingly,  the  tall  end  of  the  solution  sequence  will  be  increasing: 


ut<u.,<...<u  ( where 

J  J+1  n 

in  the  subsequence  ui»ui+i> 


1  <  j).  Furthermore,  there  are  no  sign  changes 
. ,Uj,  and  Uj  <  0.  But  it  is  possible  that 


u^  £  0.  This  means  that  a  spurious  sign  change  may  occur  in  the  tail  end  of 


the  sequence.  Of  course,  all  of  these  things  occur  in  the  Sturm  sequence, 


since  S  =  u  /t  .  Thus,  there  are  N  (A)  sign  changes  in  the  subsequence 

Sq.S^,  . . . ,S^;  there  are  no  sign  changes  in  the  subsequence  S^ , +  ^ , . . . , Sy 

S.  <  S,  ,  <  ...  <  S  ;  and  S.  <  0.  Also  S  ,  >  0,  since 
j  J+1  n’  J  n+1 

sgn  Sn+j  =  sgn  u(A),  for  small  h. 

We  now  have  two  possibilities:  either  S  <  0  or  S  >  0.  If  S  <  0, 

n  n  n 

then  a  sign  change  occurs  between  S  and  S  ..  Therefore, 

~  n  n+1 

c(h,A)  =  Nq(A)  +  1  =  N(A).  If  S^  >  0,  then  a  single  sign  change  occurs  in 
the  subsequence  Sj,  S^  +  ^ ,  .  .  .  ,  S^+^ .  Again  c(h,A)  =  Nq(A)  +  1  =  N(A).  This 


concludes  Case  2. 
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Until  now,  we  have  assumed  that  /3  (X)  >  0.  If  8„(X)  s  0,  then  we  must 

o  0 

delete  the  first  row  and  column  of  the  matrix  A(X)  in  (3.8),  and  set 
uQ  =  0.  Theorem  4.2  is  changed  slightly  as  Indicated  in  Remark  1.  After 
these  modif ications,  the  proof  is  the  same  as  above.  Q.E. D. 

Remark  3.  In  Theorem  4.3,  X  is  not  allowed  to  be  an  eigenvalue  of  (1.1). 

If  X  =  X^  is  an  eigenvalue,  then  for  small  e  >  0,  and  small  h,  the  dis¬ 
crete  problem  has  a  unique  eigenvalue  X^  in  (X^-e,X^+e).  (See  Theorem  A. 1 
in  the  appendix. )  If  we  do  not  use  quadratures  in  the  finite  element  discre¬ 
tization,  then  X^  <  X^.  If  this  were  also  true  with  quadratures,  then 
Theorem  4.3  would  be  valid  for  X  =  X^.  Unfortunately,  this  does  not  seem  to 
be  true,  in  general.  Moreover,  if  we  were  to  use  a  pure  finite  element  dis¬ 
cretization  without  quadratures,  then  the  matrix  A(X)  in  (3.8)  would  not  have 
the  structure  necessary  for  the  application  of  Sturm  sequences. 

Lemma  4. 1. 

(1)  Let  A^  <  X'  <  X"  <  A^.  If  the  interval  tX'.X"]  contains  no 
eigenvalues  of  (1.1),  then  N(X')  =  N(X"). 

(2)  If  X^  is  an  eigenvalue  of  (1.1),  then  for  small  c  >  0, 

N(Xk-c)  =  N(Xk)  and  N(Xk+e)  =  N(Xk)  +  1. 

Proof.  (1)  For  small  h,  the  discrete  problem  has  no  eigenvalues  In 
[X',X"J.  (See  the  appendix  for  properties  of  the  approximating  eigenvalues 
A  )  Theorem  4.3  implies  that  c(h,X)  =  N(X')  and  c(h,A")  =  N(X"),  for 
small  h.  Theorem  4.1  (1)  implies  that  c(h.X')  =  c(h,X"),  and  therefore 
N(X')  =  N(X") . 

(2)  We  shall  consider  two  cases,  according  to  the  possibilities  that 

u(  1 ,  X.  )  =  0  or  not. 
k 
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Case  1 .  ud.A^)  *  0. 

In  this  case,  the  number  of  zeros  of  u(x,A)  in  (0,1)  does  not  change  if  A 

is  moved  slightly  away  from  A^.  Therefore  Nq(A^-c)  =  NQ(Ak)  =  NQ(Ak+e).  By 

definition,  <r( A.  )  =  0,  so  N(A  )  =  N  (A.  ).  For  small  h,  the  discrete 
K  K  O  K 

problem  has  a  unique  eigenvalue  A^  in  the  interval  (A^-c,A^+c)  (see  the 
appendix).  Therefore  c(h,A^+e)  =  c(h, A^-c)  +  1.  Theorem  4.3  implies  that 
N(Ak+e)  =  N(A^-e)  +  1.  This  implies  that  cr(Ak+e)  =  1  and  c-fA^-c)  =  0,  so 
N(AR-e)  =  N(Ak)  and  N(Ak+e)  =  N( AR)  +  1. 

Case  2.  u(l,Ak)  =  0. 

In  this  case,  3^(A)  =  0,  <r(A)  =  0  and  N(A)  =  Nq(A),  for  A^  <  A  <  A^. 

When  A  is  moved  slightly  from  Ak  to  Ak±e,  the  number  of  zeros  of  u(x,A) 
in  (0,1)  either  remains  the  same,  or  increases  by  1.  Thus  each  of  the 


numbers  N  (A  +c),  N  (A  -c)  equals  either  N  (A  )  or  N  (A  )+l.  As  in 

OKOiC  OK  OK 

case  1,  c(h,Ak+e)  =  c(h,Ak~c)+l,  for  small  h,  and  therefore 

N  (A  +e)  =  N  (A  -e)  +  1.  This  implies  that  N  (A,-e)  =  N  (A  )  and 
ok  ok  ok  ok 

N  (A,  +e)  =  N  (A.  )  +  1.  Q.  I 

ok  ok 

Lemma  4.2. 


Q.E.D. 


(1)  If  (1.1)  satisfies  the  limit  assumption  (L3),  then  Sj(h,A+)  2  0 
for  all  h  =  1/n,  and  1  <  1  S  n  +  1. 

(2)  If  (1.1)  satisfies  the  limit  assumption  (L2),  then  there  exists 

hQ  >  0  and  AQ  in  (A^.A^),  such  that  S^Oi.Aq)  t  0  for  all  h  =  1/n  <  h^, 
and  1  <  i  <  n  +  1 . 


Proof.  (1)  In  Greenberg  (7,  Lemma  3.3]  it  Is  shown  for  a  matrix  A  of  type 
(3.8),  that  if  Z  0  for  0  S  1  S  n-1,  and  a,.  :£  0  for  0  <  1  S  n, 

then  detfA]  >  0.  Referring  to  (3.9),  we  see  that  this  implies  (1). 
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(2)  Assumption  (L2)  implies  that  lim  q  (A)  =  -oo.  Therefore  there  is  a 

A-*Aj 

* 

Aq  in  (A^.A^)  such  that  q  (A^)  <  0.  Referring  to  (3.9),  we  see  that 

a^fh.Ag)  <  0  for  all  h  =  1/n,  and  1  <  i  <  n-1,  while  b^h.A)  >  0  for 

all  A,  h  =  1/n  and  0  <  i  <  n  -  1.  We  must  still  deal  with  and  a  . 

U  n 

(Recall  that  if  0O(A)  3  0,  then  the  first  row  and  column  of  A(A)  in  (3.8) 

are  omitted,  so  a^  does  not  occur.  Similarly,  if  0^(A)  =  0,  then  a^ 

does  not  occur.  Therefore  we  may  assume  that  0n(A)  >  0  and  0.(A)  >  0.) 

2  u  i 

Define  aQ(h,A)  =  jj-q0(A),  bQ(h,A)  =  Pl(A)  -  hg|A|p0(A) , 

2  2 

a  (h.A)  =  jp-q  (A)  and  b  .  (h,A)  =  p  i ) ( A )  +  h^y^yp  (A).  Then 
n  n  n-i  n —  piia;  n 

2 

a0(h,A0)  <  0,  an(h,A0)  <  0,  bQ  -  aQ  =  bQ  -  aQ  and  =  Vl'V  We 

must  show  that  b„(h,A„)  >  0  and  b  ,(h,A„)  >  0  for  small  h.  Let 
0  0  n-1  u 

m  =  max  p(x,A  ),  and  recall  that  p(x,A)  >  k  >  0,  by  assumption  (S2).  Then 
0<x<l  ° 

Eolh-V  s  k'hl^TlJTl"  5n-i(h-xo>  *  k'hlgTT^Tlm'  Let 

ho  =  "‘"<s's^  s'tw1’-  The"  Vh’V>0  a"d  Vi(h-V>0'  for 


0  <  h  <  hQ. 


A  Second  Proof  of  the  Shooting  Theorem. 


Q.E.D. 


(1)  Suppose  that  (1.1)  satisfies  the  monotonicity  assumptions  (M1)-(M4). 

Lemma  4. 1  shows  that  N(A)  is  a  piecewise  constant  function,  with  jump  dis¬ 
continuities  at  the  eigenvalues  A^,  and  N(A^-e)  =  NtA^),  N(Ak+e)  =  N(Ak)+l. 

Therefore  N(A")-N(A')  equals  the  number  of  Jump  discontinuities  in 

[A', A"),  which  equals  the  number  of  eigenvalues  in  [A', A"). 

(2)  By  definition  of  the  eigenvalue  A^,  utx.A^)  has  exactly  k-1 
interior  zeros.  By  definition  of  N(A),  N(A^)  =  k-1.  For  A^_^  <  A  <  A^, 

N ( A )  =  N(Ak)  because  the  interval  ( A , A^ )  contains  no  eigenvalues  of  (1.1). 
Thus  N(A)  =  k  -  1  for  A^  <  A  S  A^.  Now  suppose  that 

N(A')  =  j  <  k  <  N( A" ) .  Then  the  interval  [A', A")  contains  exactly  k- j 


’4. 


eigenvalues  and  A  £  A  <  ...  <  A.  <  A 


(3)  Suppose  that  (1.1)  also  satisfies  either  limit  assumption  (L2)  or 
(L3).  By  Lemma  4.2,  there  exists  A^  such  that  Sj(h,AQ)  £  0  for  small 
h  =  1/n,  and  1  <  1  £  n+1.  This  implies  that  c(h,Ag)  =  0.  If  is  not 

an  eigenvalue  of  (1.1),  then  c(h,Ag)  =  N(AQ),  for  small  h,  and  therefore 
N(Aq)  =  0. 

Suppose  that  A^  is  an  eigenvalue  of  (1.1).  For  small  h,  the  discrete 
problem  has  a  unique  eigenvalue  A^  in  (Aq-c.Aq+g)  (or  in  [Ag,A0+c),  if 
Aq  =  Aj).  Theorem  4.1  (1)  Implies  that  c(h,AQ+e)  =  c(h, A^-e)  +  1.  But  since 
c(h,A)  is  a  nondecreasing  function  of  A,  and  since  c(h,AQ)  =  0,  it 
follows  that  c(h,AQ-c)  =  0.  Therefore  c(h,AQ+e)  =  1,  for  small  h.  Since 
c(h,AQ+c)  =  N(Aq+g),  this  implies  that  N(Aq+c)  =  1.  By  Lemma  4.1, 

N(Aq+c)  =  N(AQ)  +  1,  therefore  N(AQ)  =  0.  Now  (1)  implies  that  (1.1)  has  no 
eigenvalues  in  (A^.A^),  and  N(A)  =  0  for  Aj  <  A  S  A^.  If  A  >  AQ,  then 
any  eigenvalues  in  (A^.A)  are  contained  in  (A^.A),  and  the  number  of  these 


eigenvalues  is  N(A)-N(AQ)  =  N(A). 


Q.E.D. 


Sturmian  Theorems  for  the  Discrete  Problem. 


In  this  section,  we  shall  prove  discrete  versions  of  the  Sturm  compari¬ 
son,  oscillation  and  separation  theorems.  For  1  £  i  £  n,  let  A^(A)  be  the 
ixi  matrix  which  is  the  upper  left  corner  of  the  matrix  A ( A )  in  (3.8): 


(5.1)  A  (A)  = 


=  detfA^CA)]  is  the  itb  function  in  the  Sturm  sequence  for  A(A).  The 
Sturm  sequence  for  A^ ( A )  is  Sq.Sj, S^, . . . ,S^ . 

Definition  5. 1 .  For  1  S  1  S  n,  and  for  a  given  A  in  (A, ,  A_ ) ,  c  (A)  is 
the  number  of  sign  changes  in  the  sequence  Sq( A) , S^ ( A) , . . . , S^ ( A) ,  after  the 
zero  terms  (if  any)  have  been  omitted  from  the  sequence. 

(Note  that  c^(A)  has  been  defined  for  1  <  i  S  n.  We  shall  not  use  the 

notation  c  ,(A)  for  c(A),  or  A  ,(A)  for  A(A).)  The  following  theorem 
n+ 1  n+i 

is  proved  in  Greenberg  [7,  Theorem  3.2]. 

Theorem  5.1.  Suppose  that 

(1)  a. (A)  is  a  strictly  increasing  function,  for  0  £  j  £  1-1, 

3 

(2)  b  (A)  is  a  nonincreasing  function,  for  0  i  J  S  1-1, 

J 

(3)  bj(A)  has  no  zeros,  for  1  <  J  <  i-2. 

Let  A'  <  A"  be  numbers  in  (A^.Ag).  Then  det[A^(A)]  has  exactly 
c,(A")  -c.(A')  different  zeros  in  the  interval  [A', A"). 


( bo”ao2  3 


-bo  (Vbrai3  -bi 


-bi  ‘VW  _b2 


“bi-3  (bi-3+bi-2~al-2)  bi-2 


-b 


i-2 


(bi-2+bi-l  ai-l3 


In  the  present  situation,  where  A(A)  is  the  finite  element  matrix  corres¬ 
ponding  to  a  discretization  of  (1.1),  b.  =  p  i  >  0,  so  condition  (3)  in  the 

J  )*2 

theorem  is  satisfied.  If  (1.1)  satisfies  the  monotonicity  assumptions 
(Ml)  -  (M3),  then  conditions  (1)  and  (2)  are  also  satisfied.  Note  that 
assumption  (M4)  is  not  required,  because  otj(A)  and  fJj(A)  appear  only  in 
the  last  row  of  A(A),  so  they  do  not  appear  in  any  of  the  matrices  A^(A), 

1  <  1  <  n. 

Let  B  =  B ( A )  be  the  nx(n+l)  matrix  obtained  from  A(A)  by  deleting 
the  last  row: 


(5.2)  B  = 


T 

Let  u  =  u(A)  =  (u  (A) , u, ( A) . u  (A))  be  the  solution  of  the  discrete 

0  1  n 

version  of  the  unitial  value  problem  (1.4).  In  other  words,  u  is  the  solu¬ 
tion  of 

(5.3)  Bu  =  0,  uQ  =  |3  (A). 

We  shall  usually  assume  that  £q(A)  >  0-  with  an  occasional  remark  on  the 
case  /3Q(A)  =  0-  Although  B  and  u  depend  on  the  mesh  size  h  as  well  as 
A,  we  shall  usually  suppress  h  from  the  notation,  since  the  mesh  size  will 
usually  be  fixed,  in  this  section. 

The  last  n  -  1  equations  of  (5.3)  are 


(5.4)  — bi_iui_i  ( b  i  _  i +bi — ai )  ui  —  t)iui  +  i  =  °*  for  1  ^  i  S  n-1. 

This  implies  that  two  consecutive  terms  u^>ui+j  cannot  both  be  zero.  (We 

have  already  discussed  this  after  equation  (4.5),  which  coincides  with  (5.4).) 

Also,  (5.4)  implies  that  if  u.  =  0,  then  u.  . u.  ,  <  0  (since  all  b.  >  0). 

i  i-1  i+1  j 


Definition  5. 2.  We  shall  say  that  u(A)  has  a  sign  change  at  u^(A)  If 

either  Uj_j(A)Uj(A)  <  0,  or  Uj_^(A)  =  0. 

Note  that  if  u  =  0,  then  j-1  >  0  and  u  u  <  0. 

J  *  J  ^  J 


Theorem  5.2  (First  Comparison  Theorem).  Suppose  that  (1.4)  satisfies  the 
monotonicity  conditions  (Ml)  -  (M3).  Let  Aj  <  A^  be  numbers  in  (A^.A^). 
Then: 

(1)  u(A2)  has  at  *east  as  many  sign  changes  as  u(Aj).  Furthermore, 
u(A^)  has  more  sign  changes  than  u(A^)  if  and  only  if  u^(A)  has  a  zero  in 
the  interval  (A^.A^)- 

(2)  If  the  ith  sign  change  for  u(A^)  occurs  at  u^  ( A ^ ) ,  and  for 

u(A„)  at  u.  ( A„ ) ,  then  j  S  J  .  Furthermore,  J_  <  j.  if  and  only  if 
2  J2  2  2  1  2  l 

u,  (A)  has  a  zero  in  (A  ,A  ). 

Ji  12 


Proof. 

(1)  By  Theorem  4.2,  u(A)  has  the  same  number  of  sign  changes  as  the 

Sturm  sequence  S„(A),  S, (A ),...,  S  (A)  for  A  (A).  Therefore  u(A)  has 
0  1  n  n 

exactly  c^lA)  sign  changes.  By  Theorem  5.1,  c^(Aj)  <  c^tA^)-  Furthermore, 

c  (A.)  <  c  (A_)  ix  and  only  if  S  (A)  =  u  (A)/t  (A)  has  a  zero  in  (A  A  ). 
nln2  nnn  12 

(2)  Similarly,  the  sequence  Ug(A),  u^ (A) , . . . , Uj ( A )  has  exactly  c.(A) 

sign  changes.  By  assumption,  c.  (A.)  =  i.  Since  c.  (A.)  <  c.  (A  ),  the 

Ji  1  Ji  i  Ji  2 

sequence  u_(A„) , u. (A_) , . . . , u  ,  (A„)  has  at  least  i  sign  changes.  Therefore 
0  2  1  2  Ji  2 


J2  £  Jj.  Furthermore,  J2  <  J1  If  and  only  if  uQ(X2),  U1  ( ^-2 ) . Uj  (X,,) 

has  more  than  i  sign  changes.  The  latter  is  true  if  and  only  if  u,  (X) 

Ji 

has  a  zero  in  [X,,X_).  However,  since  a  sign  change  occurs  at  u,  (X.), 
i  &  Ji  1 

Definition  5.2  Implies  that  u.  (X.)  *  0.  Therefore  the  zero  must  occur  In 

Ji  1 


(W 


Q.  E.  D. 


Remark.  In  Greenberg  [7],  Sturm  sequences  are  discussed  from  an  axiomatic 
standpoint.  It  is  noted  there  that  one  of  the  axioms  implies  the  following 
fact:  If  S  (X  )  =  0,  then  S  (X)Si+1(X)  >0  for  XQ  -  c  <  X  <  XQ,  and 

Si(X)Si  +  1(X)  <  0  for  Xq  <  X  <  Xq  +  e.  Thus,  as  X  increases  past  XQ,  a 
sign  change  is  generated  between  S^(X)  and  S1+^(X).  Of  course,  the  same  is 
true  for  u^lX)  and  ui+1(X).  If  u^(X)u1+^(X)  <  0,  then  as  X  increases, 
Ui+1(X)  cannot  change  sign  before  u^(X)  does.  This  means  that,  as  X 
increases,  the  sign  changes  move  from  right  to  left.  Also,  if 
Ui-1  (X )ui  (X)  >  0  and  u^CXju^CX)  >  0,  then  u^X)  cannot  have  a  sign 
change  before  both  u^lX)  and  ui+1(X).  F°r  ***  Uj(Xq)  =  °*  then 
Ui-1  ^0^Ui  +  l  ^0^  Thls  means  that  sign  changes  cannot  appear  spontaneous¬ 

ly  in  the  middle  of  the  sequence.  They  must  start  at  u^iX),  and  move  to  the 
left.  This  is  exactly  what  happens  to  the  zeros  of  the  solution  u(x,X)  of 
the  continuous  problem  (1.4). 

The  following  lemma  is  proved  in  Greenberg  [7,  Lemma  3.2).  (There  is  a 
slight  difference  in  notation,  for  two  reasons.  First  of  all,  the  first  row 
of  A ( X )  is  here  indexed  by  1=0,  while  in  [7],  it  is  indexed  by  1=1. 
Second,  we  assume  here  that  none  of  the  b^  have  zeros,  while  in  [7],  this 
assumption  is  made  only  for  the  first  n-1  values  bj. ) 


Lemma  5.1.  Suppose  that,  for  0  <  j  <  n  -  1, 


(1)  a  (X)  is  a  strictly  increasing  function, 
<J 
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(2)  bj(A)  *s  a  nonincreasing  function, 

(3)  bj(A)  has  no  zeros. 

2  S 

Then,  for  1  £  1  <  n,  b,  ,  -  b,  ,  is  a  strictly  decreasing  function  on 

1  —  1  1  “  1 

any  interval  which  contains  no  zeros  of  S^(A). 

As  before,  condition  (3)  is  automatically  satisfied  here,  and  conditions 
(1)  and  (2)  are  implied  by  the  monotonicity  assumptions  (Ml)  -  (M3).  We  shall 
use  the  following  notation: 

(5.5)  u'  =  (ui-uil)/h,  for  1  S  i  £  n. 

Theorem  5,3  (Second  Comparison  Theorem).  Suppose  that  (1.4)  satisfies  the 

u' (A) 

monotonicity  conditions  (Ml)  -  (M3).  Then,  for  1  5  i  <  n,  p  i(A) — pr-y  is  a 

>-i  uiU) 

strictly  decreasing  function  on  any  interval  which  contains  no  zeros  of 
u. (A) . 


Proof.  By  Theorem  4.2,  u^  =  t^,  where  t1  =  Po/^i-lt>i-2‘ ’  ‘  ^0^ '  There- 

Si  _i  . 

i-i-sr  and 


fore  Ul  —  =  b 


Ul 


p  i (A)—  =  Ul~1l  =  [" i  -  likilj  =  ifh  -b  2  ^—1 .  The  theorem 

i ui  h  [  ui  J  h  [  Ui  J  h ^  i-1  i-1  S,  J 


now  follows  from  Lemma  5.1. 


Q.  E.D. 


Remark.  If  /3q(A)  2  0.  then  the  first  row  and  column  of  A(A)  in  (3.8)  must 

be  deleted,  and  uQ  =  0.  The  conclusion  In  Lemma  5. 1  is  changed  to: 

2  S 

b.  -  bf  is  a  strictly  decreasing  function.  (This  now  agrees  with  the 

1  1  Di  / 

u; 


notation  in  [7]. )  In  this  case,  b  -  b 


2Si_,  _  ,  “i  +  1 


1  i  S 


=  hp 


1+2Ui  +  l 


,  so  the  conclu¬ 


sion  in  Theorem  5.3  remains  the  same,  except  that  2  i  i  ^  n.  (For  i  =  1, 

uj  b 

Pi —  =  which  is  a  nonincreasing  function.) 

-  u,  h 
2  1 


The  following  lemma  is  implied  by  Lemmas  3.4  and  4.2  in  Greenberg  [7], 


Lemma  5.2.  Suppose  that  (1.1)  satisfies  assumptions  (M2)  and  (LI).  Then 
11m  S  (A)  =  (-1 )  ioo,  for  l$li  n+1. 

"V  .  A  * 


Theorem  5. 4  (Oscillation  Theorem). 

(1)  Suppose  that  (1.1)  satisfies  (Ml)  -  (M4)  and  (LI).  Then  there  is  a 
number  h^  >  0  and  am  Integer  m,  such  that  for  mesh  size  h  <  hg,  and 


h  =  1/n,  the  discrete  problem  has  eigenvalues  A  <  A  , 

m  m+1 


<  ...  <  A 


(2)  Suppose  that  (1.1)  satisfies  (Ml)  -  (M4),  (LI)  and  either  (L2)  or 
(L3 ) .  Then,  for  any  mesh  size  h  =  1/n,  the  discrete  problem  has  eigenvalues 

xi  <  *2  <  <  Vi¬ 

no 

In  (1)  and  (2),  the  eigenvector  u  ,  corresponding  to  A^,  has  exactly 
k  -  1  sign  changes. 


Proof. 

(1)  If  u(A)  is  a  solution  of  the  discrete  initial  value  problem  (5.3), 
then  u(A)  has  exactly  cn(A)  sign  changes.  If  A  is  an  eigenvalue  of 
(1.1),  then  u(A)  also  satisfies  the  last  equation  in  (3.7),  and 
det [ A(A) ]  =  0. 


Let  Aj  <  Aq  <  A^,  and  suppose  that  A^  is  not  an  eigenvalue  of  (1.1) 
and  N(Aq)  =  m-1.  Theorem  4.3  implies  that  there  is  a  number  hQ  >  0.  such 
that  when  h  <  hQ,  then  c(h,Ag)  =  N(AQ)  =  m-  1.  We  may  suppose  that  h^ 
is  small  enough,  so  that  h^  <  1/m.  In  the  following,  we  assume  that  h  <  h^. 
By  Lemma  5.2,  11m  S^(h,A)  =  ( —  1 ) 1  oo.  This  implies  that  for  A  near  A^, 

c(h,A)  =  n+1.  Since  c(h,Ag)  =  m-1.  Theorem  4.1  (1)  implies  that  the 
discrete  problem  has  c(h, A)  -  c(h, A^)  =  n-m  +  2  eigenvalues  in  the  inter¬ 
val  [A-,A).  We  shall  denote  them  by  A  <  A  ,<...<  A  ,,  where  A_<  A 
0  •'mm+l  n+1  0m 

and  A  ,  <  A.  Since  A,  is  the  only  eigenvalue  in  [A.  ,A.  1,), 
n+1  k  jo  k  k+1 


c(h,X^+j)  =  cfh.A^)  +  1.  Also  c(h,An+j)  =  c(h,A)  -  1  =  n.  This  implies  that 

c(h,A^)  =  k-1,  for  m  i  k  <  n+1.  At  any  eigenvalue  A^, 

Sn+^(h,A^)  =  det[A(h, A^) ]  =  0.  Therefore  cCh.A^)  =  c^fh.A^).  This  shows 

(k) 

that  the  eigenvector  u  =  u(A^)  has  exactly  k-1  sign  changes. 

(2)  If  (1.1)  also  satisfies  (L2)  or  (L3)  then  Lemma  4.2  implies  that  A 
can  be  chosen  so  that  c(h,Ag)  =  0,  for  small  h.  In  this  case,  m  =  1. 

Q.E.D. 

We  shall  conclude  this  section  with  a  discussion  of  the  Sturm  separation 
theorem.  This  theorem  does  not  involve  a  parameter  A.  It  asserts  that  if 
u(x)  and  v(x)  are  linearly  Independent  solutions  of 

(5.6)  (p(x)y')'  +q(x)y  =  0, 

then  the  zeros  of  u(x)  and  v(x)  are  interlaced.  Let  C  be  the 
( n- 1 ) x ( n+ 1 )  matrix,  corresponding  to  A( A)  in  (3.8),  with  the  first  and 


Cn-2  "bn-2 

bn-2  Cn-1  bn-l 

where  c^  =  b^  ^  +  b^-a^.  Here,  the  coefficients  b^.c^  are  constants, 
with  b^  >  0.  The  discrete  problem  corresponding  to  (5.6)  is 

(5.8)  Cy  =  0, 

T 

where  y  =  (yQ, yj * ■ • • • yn)  •  As  before,  if  y  is  a  nontrivial  solution  of 

(5.8) ,  then  two  consecutive  terms  y^.y^  cannot  both  be  zero.  If  y^  =  0 
(where  1  <  1  <  n-1),  then  yj-^yj  +  j  <  However,  in  the  present  sltua- 


last  rows  omitted: 


(5.7)  C  = 


-bo  C1 

-b. 


C2  "b2 


-b 


n-3 


3! 


tlon,  the  possibility  that  y^  =  0  is  not  excluded.  We  continue  to  use 
Definition  5.2  for  the  position  of  a  sign  change.  In  particular,  if  y^  =  0, 


then  y  has  a  sign  change  at  y^. 


Definition  5.3.  The  discrete  Wronskian  of  u  =  (u_.u . u  )  and 

y  j  n 

T 

V  =  (V  V . v  )  is  w(u,v)  =  (w  ,w. . w  ),  where 

u  l  n  o  l  n-l 

wi  =  uivi+i "  viui+r  for  0  *  1  *  n'1- 


T  T 

Lemma  5.3.  Let  u  =  (u„.u . u  )  and  v  =  (v_.v . v  )  be  solutions 

-  0  1  n  0  1  n 

of  (5.8),  and  let  w(u,v)  =  (w^.w . w  .). 

0  1  n-1 


(1) 

u 

and 

v  are  linearly  dependent  if  and  only 

if 

w(u,v)  =  0. 

(2) 

bi 

Wi  = 

Vo* 

for  1  £  1  <  n-  1. 

(3) 

If 

u 

and  v 

are  linearly  Independent,  then 

W1 

*■  0  for 

0  <  i  <  n 

-  1 

,  and  the 

w^  all  have  the  same  sign. 

Proof. 

(1) 

We 

may 

assume 

that  neither  u  nor  v  is  the 

zero 

vector.  If 

and  v  are  linearly  dependent,  then  there  Is  a  constant  c  such  that 
v  =  cu.  Then  wJ  =  u^,  v^  +  1  -  v^u^  =  u^cuJ  +  1  -  cUjUj  + ^  =  0.  Conversely,  sup- 


pose 

that 

w(u,v)  = 

0 .  Then 

UiVi+l " ViUi+l 

0, 

for  0  S 

i  S  n  - 

1.  If 

v,  * 

0  for 

0  <  i  < 

n,  then 

Ul  =  Ul+1.  Thus 

u0 

-  Hi  = 

=  Un  _ 

c,  and 

1 

Vl  vUi 

v0 

Vl 

vn 

u  =  cv.  If 

Vi  =  °’ 

then  the 

relation  uiv^+^ 

~Vi 

u.  =  0 

l  +  l 

impl ies 

that 

u,  =  0,  since  v,  ,  *  0.  (Or,  if  1  =  n,  u  ,v  -v  .u  =0  implies  that 
i  i+1  n- Inn- In 

u  -  0.  )  Let  c  =  (or  if  1  =  n,  c  =  ^S^),  and  let  y  =  u  -  cv. 

n  vul  vn.j 

Then  y  is  a  solution  of  (5.8)  and  two  consecutive  terms  y^.y^  are  zero- 
Therefore  y  =  0  and  u  =  cv. 

(2)  The  equations  (5.8)  for  u  and  v  are  of  the  form 
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£v. 


v  V 


A lViVvvvvv V  •C' 


•  „  ^  . 


v:.v 


"Jl  ">  ’.K  ".X  "J  '_K  TV'  IT* 


(5.9) 


‘bJ-lUJ-l  *  Yj  -  bjuj»l  ■  °' 
;bj-ivj-i  *  cjvj  ■  bjvj*i  *  °- 


Multiply  the  first  equation  in  (5.9)  by  v,, 

J 

and  then  subtract.  This  gives  the  equation  b 

ttot  Yj  '  bj-l"j-l  ‘  bJ“2WJ-2  '  '  • '  =  b0K0- 
(3)  This  follows  immediately  from  (2). 


and  the  second  equation  by  u^, 
*  ^JWJ’  which  implies 


Q.E.D. 


T 

Theorem  5. 5  (Separation  Theorem).  Let  u  =  (u„,u, ,....u  )  and 

0  1  n 

T 

v  =  (Vq.Vj, . . . , v  )  be  linearly  Independent  solutions  of  (5.8).  Suppose  that 
u  has  a  sign  change  at  u^  and  at  u^  (where  i  <  J),  and  no  sign  changes 
between  these.  Then  v  has  a  sign  change  at  some  v  ,  where  i  <,  k  S  j. 


Proof.  Since  u  cannot  have  a  sign  change  at  uQ,  0  <  i.  The  sign  changes 
at  Uj  and  u^  imply  that  u^jUj  <0,  ui  *  0,  uj_iuj  5  0  and  uj  *  °- 
Since  u  has  no  sign  changes  between  u^  and  u^,  we  may  suppose  that  the 


terms  Uj.u  ,  . 

.  ••uj_2  are  all  positive  u^  <  0,  u^  ^  > 

0  and 

UJ  < 

Suppose  that  v 

has  no  sign  change  in  the  Indicated  range. 

Then  we 

may 

suppose  that  v^ 

_j,v^ . VJ-2’Vj-l  are  a1^  positive,  and 

6 

Al 

> 

Let 

w(u,v)  =  ( wq »  w 1 , 

••"Vi1-  Then  ui-i  =  ui-ivi  -  vi-iui  <  0 

and 

w.  . 


u,  ,v  -  v  u,  >  0.  This  contradicts  Lemma  5.3  (3). 


Q.E.D. 


■VLVWUvysvvvx  vv -•> .< v***' ■>  •> 'p-prp  v -Jr?" "f.-p.-p.vjv 


6.  An  Incorrect  Algorithm. 


We  shall  consider  an  alternative  algorithm  for  calculating  the  n 
eigenvalue  of  (1.1).  Let  u(x,A)  denote  a  solution  of  the  initial  value 
problem  (1.4),  and  let 


(6.  1) 


u(x,A)  =  otj(A)u(x,A)  +  0^ (A)u' (x, A) . 


It  might  seem  natural  to  count  the  zeros  of  u(x,A)  (for  fixed  A),  rather 
than  counting  the  zeros  of  u(x,A),  with  a  correction  using  u(l,A),  as  in 
§2.  In  other  words,  if 

(6.2)  N'(A)  equals  the  number  of  zeros  of  u(x,A)  in  the  interval  (0,1), 
we  might  conjecture  the  following: 


Hypothetical  Theorem.  If  (1.1)  satisfies  the  monotonicity  assumptions 

(Ml)  -  (M4) ,  and  either  the  limit  assumption  (L2)  or  (L3),  then  for 

A^  <  A  <  A^,  (1.1)  has  exactly  N'(A)  eigenvalues  in  the  interval  [A^.A). 


This  would  lead  to  an  alternative  shooting  method  for  finding  A^.  However, 
the  above  Hypothetical  Theorem  turns  out  to  be  false.  Indeed,  it  is  already 
false  for  linear  eigenvalue  problems.  It  depends  on  monotonicity  properties 
of  eigenvalues,  which  are  valid  for  Dirichlet  boundary  conditions,  but  not  for 
general  boundary  conditions.  Because  of  this,  the  alternative  shooting  method 
does  not  always  work.  We  shall  presently  given  an  example  where  the  alterna¬ 
tive  shooting  method  fails.  This  section  has  been  included,  because  we 
believe  that  it  is  useful  to  point  out  that  a  numerical  method  can  fail.  It 
is  especially  important  for  this  method,  which  seems  to  attract  believers 
easily,  and  which  may  be  used  in  applications. 

Note  that  if  u(xQ,A)  =  0,  then  A  is  an  eigenvalue  on  the  interval 


SO 


a 
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[O.Xq],  and  xQ  is  a  critical  length  (as  used  in  the  invariant  Imbedding 
method).  Thus  N' (X)  is  the  number  of  critical  lengths  which  correspond  to 
A  (and  which  are  less  than  1). 

We  shall  now  indicate  the  role  played  by  monotonicity  properties  of 
eigenvalues.  We  shall  consider  linear  eigenvalue  problems.  In  this  case, 
(1.1)  has  the  form 


(6.3) 


’(p(x)u' )'  +  (Ar(x)-q(x) )u  =  0,  for  0  5  x  5  1, 
aQu(0)  +  0ou' (0)  =0, 
a^u(  1 )  +  /SjU'd)  =  0, 


where  r(x)  >  m  >  0,  and  0^,0^  are  constants  (for  1  =0,1).  For 
0  <  y  <1,  let  An(y)  denote  the  n^  eigenvalue  of  the  problem: 


(6.4) 


(p(x)u')'  +  (Ar(x)-q(x) )u  =  0,  for  0  5  x  5  y, 
aQu(0)  +  0Qu' (0)  =0, 
a^u(y)  +  0^u' (y)  =  0. 


For  given  n,  we  ask  the  following  question: 


(Q)  Is  An(y)  a  decreasing  function  on  (0,1)? 


If  the  answer  to  (Q)  is  yes,  for  all  n,  then  the  Hypothetical  Theorem  is 
true  and  the  alternative  algorithm  is  correct.  To  verify  this,  let  u(x,Aq) 
be  the  solution  of  the  initial  value  problem 


(6.5) 


(p(x)u' )'  +  (AQr(x)-q(x) )u  =  0,  for  05x51, 
u(0)  =  0O>  u'  (0)  =  -0Cq, 


and  let  u(x,Aq)  =  ctjUtx,  AQ)  +  0^u'  (x,  AQ) .  If  x^  is  a  zero  of  u(x,AQ)  in 
(0,1),  then  Aq  is  an  eigenvalue  on  [0,Xq].  In  other  words,  A^  =  ^(Xq^ 


for  some  k.  Since  A^(y)  is  a  decreasing  function, 

Ak  =  A^d)  <  Ak(xQ)  =  *q*  In  this  way,  each  zero  of  u(x,Aq)  corresponds  to 

an  eigenvalue  less  than  Aq.  This  shows  that  the  Hypothetical  Theorem  and 

alternative  algorithm  are  correct,  if  A^(y)  is  a  decreasing  function,  for 

all  n.  Unfortunately,  this  is  not  true  for  general  boundary  conditions. 

If  the  boundary  condition  at  the  right  endpoint  is  a  Dirichlet  condition: 

u(y)  =  0,  then  the  classical  monotonicity  theorem  tells  us  that  A^ty)  is  a 

decreasing  function,  for  all  n.  But  Greenberg  [8]  has  shown  that  if  the 

boundary  condition  at  the  right  endpoint  is  not  a  Dirichlet  condition  (i.e., 

*  0),  then  for  given  >  1,  there  exist  coefficient  functions  p(x), 

q(x),  r(x)  and  a  subinterval  [a,b]  c  (0,1],  so  that  the  eigenvalues  A^(y), 

A„(y), . . • , A  (y)  are  increasing  functions  in  [a,b].  (On  the  other  hand,  for 
d  no 

given  p(x),  q(x),  r(x),  aQ,  a  ,  0^,  there  exists  nj  2:  1,  so  that  for 
n  >  n^,  An(y)  is  a  decreasing  function  on  (0,1].)  THus,  we  cannot  expect 
the  Hypothetical  Theorem  and  alternative  algorithm  to  be  correct  for  general 
boundary  conditions.  We  now  give  a  concrete  example  where  they  fail. 


Example.  For  0  <  y  <  1,  consider  the  eigenvalue  problem 

f(p(x)u')'  +  Au  =  0,  for  0  <  x  £  y, 


(6.6) 


The  energy  norm  is 


•  u(0)  +  u' (0)  =  0, 
u' (y)  =  0. 

given  by 


(6.7) 


B(  v,  v) 


-p(0)v(0)2  + 


p(x)v' (x)2dx, 


0 


and 

(6.8) 


Ax(y) 


inf 

veC1 [0, y] 


B( v, v) 

2  ’ 

iivir 
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where  llvll 


=  r'vfx 

J  r* 


)  dx.  Putting  v(x)  si,  we  find  that 


A  (y)  <  glV’V)  =  -PJ-QI, 


(6.9) 


Aj(y)  <  0. 


We  now  consider  the  two  algorithms  (given  by  Theorem  2.1  and  the  Hypothetical 
Theorem)  for  finding  the  number  of  eigenvalues  A  <  0  (for  the  interval 
[0,y]).  We  must  solve  the  initial  value  problem: 


(6. 10) 


'-(p(x)u'  )'  =  0 
u(0)  =  1,  u' (0)  =  -1. 


Denote  the  solution  by  u(x),  and  let  u(x)  =  u' (x).  We  obtain: 
-pu'  =  constant  =  p(0),  so  that 


(6.11) 


’G‘*>  -w 
u(x) ' 1 '  £  fnV 


Since  u' (x)  <  0,  u(x)  is  a  decreasing  function  (with  u(0)  =  1).  For  a 
given  y  (0  <  y  <  1),  u'(y)  <  0  and  either 
(A)  u(y)  >  0,  or  (B)  u(y)  <  0. 

In  case  (A),  NQ  =  [number  of  zeros  of  u(x)  in  (0,y)]  =0,  <r  =  1,  and 
N  =  Nq  +  «t  =  1.  In  case  (B),  Nq  =  1,  <r  =  0,  and  N  =  NQ  +  <r  =  1.  Thus  we 
see  that  the  algorithm  of  Theorem  2. 1  counts  one  negative  eigenvalue  on 
[0,y],  for  all  y  in  the  interval  (0,1). 

On  the  other  hand,  the  alternative  algorithm,  based  on  the  Hypothetical 
Theorem  counts  the  zeros  of  u(x)  =  u'(x)  in  (0,y).  Since  u(x)  <  0, 


N'  -  0,  predicting  no  negative  eigenvalues!  Here  we  have  an  example  where 


the  Hypothetical  Theorem  and  alternative  algorithm  are  incorrect. 


Remark  1 .  Greenberg  [8]  has  shown  that  In  the  above  example,  A^(y)  is  an 
increasing  function  on  (0,1].  Thus  we  have  "reverse  monotonicity"  in  this 
example! 


Remark  2.  An  error  of  only  1  in  the  integer  function  N' (A)  can  have  a 
fatal  effect,  when  the  bisection  method  uses  N'(A)  instead  of  N(A),  as 
described  in  §2.  However,  the  error  can  be  much  larger  than  1,  if  some 
eigenvalue  A^y)  oscillates  about  a  value  A  =  AQ. 

Remark  3.  The  eigenvalues  and  critical  lengths  can  be  understood  geometrical 
ly  by  introducing  polar  coordinates  in  the  phase  plane  (see  Scott,  Shampine 
and  Wing  [16]).  This  is  usually  called  the  Prufer  transformation.  Let 


(6. 12) 


u  =  u(x,A)  =  r(x,A)  cos(x,A), 
v  =  p(x, A)u' (x, A)  =  r(x,A)  sin  0(x,A). 


The  point  U(x,A)  =  (u(x, A) , p(x, A)u' (x, A) )  moves  along  a  curve  in  the 
(u.v)-plane  as  x  varies,  with  A  fixed,  or  as  A  varies,  with  x  fixed. 
Denoting  0' (x,A)  =  —  and  r'(x,A)  =  (which  conforms  to  our 

notation  u' (x,A)  =  duCx,  \)  j ,  we  ^ave 


(6. 13) 


’u'  =  r'cos  0  -  r0'  sin  0, 
(pu')'  =  r'  sin  0  +  r0'  cos  0. 


Setting  pu'  =  r  sin  0  and  (pu' )'  +qu  =  0,  equations  (6.13)  imply 


(6. 14) 


pr'  cos  0  -  pr0'  sin  0  =  r  sin  0, 
r'  sin  0  +  r0'  cos  0  =  -qr  cos  0. 


We  can  solve  for  0'  in  (6.14),  to  obtain 

(6.15)  0'  =  -[q  cos20  + -sin20]. 

P 


Also,  (6. 12)  impl ies 


(G. 16) 


taxi  8(x, A)  =  p(x,  A) 


u' (x, A) 


Equation  (E. 16)  and  Sturm’s  second  comparison  theorem  show  that  as  A 
increases,  with  x  fixed,  the  point  U(x,A)  moves  in  a  clockwise  direction. 
On  the  other  hand,  (6.15)  shows  that  as  x  increases  from  0  to  1,  with  A 
fixed,  U(x,A)  may  move  either  in  the  clockwise  or  counterclockwise  direction 
(if  q  <  0).  However,  when  U(x,A)  crosses  the  v-axis,  it  moves  in  a  clock¬ 
wise  direction. 

We  shall  now  consider  the  boundary  conditions  in  (1.1).  Let 


(6. i7 ) 


tan  Tjg ( A )  =  -p(0,  A) 


tan  7)  (A)  =  -p(  1,  A) 


<xo(A) 

^TaT’ 

ax(A) 

prnvy- 


Let  (A)  be  the  line  v  =  [tan  t)j(A)]u  in  the  (u, v)-plane  (for 
i  =  0,1).  The  boundary  condition  aQ(A)u(0)  + £Q(A)u' (0)  =0  in  (  1)  is 
equivalent  to  tan  0(0, A)  =  tan  t}q(A),  which  means  that  the  point  U(0,A) 
lies  on  the  line  HQ(A).  Similarly,  the  boundary  condition 

a  (A)u( 1)  + (A)u' ( 1 )  =  0  means  that  U(1,A)  lies  on  Hj(A).  Since  u(x,A) 
satisfies  the  boundary  condition  at  x  =  0,  we  see  that  A  is  an  eigenvalue 
of  (1.1)  if  and  only  if  U(1,A)  lies  on  H^(A).  On  the  other  hand,  x^  is  a 
critical  length  for  A  if  uix^.A)  =  (A)u(Xq, A)  +  3^ (A)u' (x^, A)  =  0,  which 
means  that  U(Xq,A)  lies  on  H^(A). 

Consider  the  trajectory  of  the  point  U(1,A),  as  A  increases  from 
to  Aq.  U( 1 , A )  always  travels  in  the  clockwise  direction,  while  H^(A) 
rotates  in  the  counterclockwise  direction,  because  of  the  monotonicity  assump¬ 


tion  (M4).  The  number  of  eigenvalues  less  than  A^  equals  the  number  of 
times  U ( 1 , A )  crosses  H^(A).  A  zero  of  u(l,A)  corresponds  to  a  point 
where  U(1,A)  crosses  the  v-axis.  By  Sturm’s  first  comparison  theorem,  such 


zeros  travel  to  the  left  in  (0,1),  and  become  Interior  zeros  of  u(x,A). 
Since  U(1,A)  never  reverses  direction,  exactly  one  eigenvalue  occurs  between 
any  two  consecutive  zeros  of  u(l,A).  The  correction  term  <r(A)  tells  us  if 
U(  1 , A )  has  crossed  H^(A)  after  the  last  crossing  of  the  v-axis.  In 
effect,  a  proof  of  the  shooting  theorem  can  be  given  in  this  framework. 

Now  consider  the  trajectory  of  U(x,Aq),  as  x  increases  from  0  to  1 
for  fixed  A  =  A^.  A  critical  length  occurs  each  time  that  U(x,Aq)  crosses 
HjfAg).  Thus  we  can  count  the  critical  lengths  by  counting  the  number  of 
times  that  tan  8(x,Aq)  =  tan  t)^(Aq).  (See  equation  (5.3)  and  inequality 
(5.4)  in  Scott,  Shampine  and  Wing  [16].)  A  zero  of  u(x,AQ)  corresponds  to  a 
point  where  U(x.A^)  crosses  the  v-axis.  As  indicated  above,  such  a  cros¬ 
sing  must  occur  in  the  clockwise  direction.  Between  two  crossings  of  the 
v-axis,  U(x,Aq)  must  cross  the  line  H^Aq).  Thus,  there  is  a  critical 
length  between  any  two  zeros  of  u(x,AQ).  However,  U(x,AQ)  may  reverse  its 
direction  several  times  between  two  crossings  of  the  v-axis.  Therefore, 
there  may  be  many  critical  lengths  between  two  zeros.  On  the  other  hand,  if 
Aj  <  Aq  <  then  U(x,Aq)  never  crosses  the  v-axis.  In  this  case,  it  may 

never  cross  H^fAp),  even  though  A^  <  A^.  This  is  what  happens  in  the 
example  (6.6).  This  explains  why  the  number  of  critical  lengths  N' (A^)  may 
be  slightly  smaller  than  N(A  ),  or  a  great  deal  larger  than  N(A  ). 


Appendix 


In  this  appendix,  we  shall  prove  some  facts  about  the  convergence  of 
approximate  eigenvalues.  The  main  facts  are  that,  as  h — >0,  the  approximate 

eigenvalues  converge  to  the  true  eigenvalues,  and  to  nothing  else.  Theorems 
A. 1  and  A. 2  make  this  statement  precise.  These  facts  were  used  in  §4,  especi¬ 
ally  in  the  proof  of  Lemma  4.1,  and  the  second  proof  of  the  shooting  theorem. 
We  are  including  these  proofs  here  for  the  sake  of  completeness.  For  a  fixed 
value  of  A,  consider  the  linear  eigenvalue  problem 

"(p(x,A)u' )'  +  q(x,A)u  +  pu  =  0,  for  0  <  x  <  1, 

(A.  1)  ■  a0(A)u(0)  +  j3  (A)u'(0)  =  0, 

a  U)u(l)  +  p^Aiu'  (1)  =  0. 

Here  the  eigenvalue  is  p.  The  weak  form  of  the  equation  is 
(A. 2)  B ( A ; u , v )  =  p<u,v>,  for  all  v  e  H*[0,1], 

r1 

where  B(A;u,v)  is  the  energy  inner  product  (3.1),  and  <u,v>  =  uvdx.  The 

J0 

monotonicity  conditions  (Ml)  -  (M4)  imply  that  B(A;u,u)  is  a  strictly 
decreasing  function  of  A.  By  the  variational  characterization  of  eigen¬ 
values,  the  k^  eigenvalue  p^(A)  is  a  continuous,  strictly  decreasing 
function  of  A.  The  k*'*1  eigenvalue  A^  of  (1.1)  is  the  unique  zero  of 
HkU). 

We  shall  distinguish  between  two  discrete  problems:  The  pure  finite 
element  discretization,  where  the  integrals  are  not  replaced  by  quadratures; 
and  our  discretization  (3.3)  and  (3.7),  where  the  integrals  are  replaced  by 
quadratures.  Each  of  these  corresponds  to  a  linear  eigenvalue  problem.  The 
weak  formulation  of  the  pure  finite  element  discretization  is:  Find  u  e  S^, 
such  that 
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(A. 3) 


B( A; u, v)  =  p,<u,v>,  for  all  v  €  S,  . 

h  h 

The  corresponding  equation  for  our  discretization  is 

(A. 4)  B.(A;u,v)  =  p  <u,v>  for  all  v  e  S,  , 

n  h  h  h 


where  Bh(A;u,v)  is  obtained  from  B(A;u,v)  by  replacing  the  integrals  by 
quadratures  (as  indicated  in  §3),  and  <u,v>  is  the  trapezoid  rule  quadra- 

r1 

ture  for  uv  dx.  If  (1.1)  satisfies  the  monotonicity  assumptions  (M1)-(M4), 
J0 

then  the  eigenvalues  p^  ^(A)  of  (A. 3)  and  p^  ^(A)  of  (A. 4)  are  continuous, 
strictly  decreasing  functions  of  A.  The  k^  eigenvalue  A^  ^  of  our  dis¬ 
cretization  (3.7)  is  the  unique  zero  of  p^  ^(A). 

The  next  two  lemmas  estimate  the  errors  Ip^AJ-p^  ^(A)|  and 


^.hU] 


k,  h 


(A)  |  . 


Lemma  A. 1.  Let  A^  <  A'  <  A"  <  A^.  For  each  integer  k  >  0,  there  is  a 
constant  C  =  C(k)  >  0,  so  that  Ip^tA)  -  p^  ^(A) I  S  Ch2,  for  A'  £  A  <  A". 

Proof.  Let  u^  =  u^(x,A)  denote  an  eigenfunction  of  (A. 2)  corresponding  to 

the  eigenvalue  p  (A)  (such  that  ||u  ||  =1).  It  is  known  that 

K  H1 (0, 1  ] 

(A. 5)  Ip  (A)-p  (A)|  <C  inf  ||u.  -v||2 

1  veS  H  [0,  1] 

h 

for  A'  <  A  <  A",  where  is  a  constant  which  depends  only  on  the  coeffi¬ 

cient  functions  in  (1.1)  (see  Babuska,  Osborn  [3]).  It  is  also  known  from 
finite  element  approximation  theory  (see  Ciarlet  [6])  that 

(A.  6)  inf  II u  -  v||  <  C  h||u  || 

veS  *  h  [0, 1 )  ^  K  h  [  0 ,  1  ] 

n 

Furthermore,  (A.l)  implies  that  =  --[  p' u'  +  (q+M,)u  ] ,  therefore 

|C  p  K  K  K 
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||u"||  <  C  for  A  <  A  £  A",  and 

L  l 0 , 1|  J 

(A.  7)  I! u  ||  <  C  for  A'  <  A  <  A". 

K  H  [0, 1 ] 

Inequalities  IA.5),  (A. 6),  (A. 7)  imply 

(A.  8)  |pk  h(A)  -  /i  (A)  |  <  Ch2,  for  A'  £  A  <  A". 

Q.E.  D. 

We  shall  now  consider  the  error  |pk  ^(A)  ~  ^(A)|,  caused  by  the 

quadratures.  We  shall  use  the  modulus  of  continuity  m^(S)  of  a  continuous 

function  f(x).  Recall  that  m  (5)  =  sup  |f(x)-f(y)|.  If  f(x)  is  con- 

|x-y|<3 

tinuous  on  a  compact  interval,  then  lim  m  (S)  =  0.  For  a  function  such  as 

SO  * 

q(x,A),  which  is  continuous  on  [0, 1 ]x[A' ,  A"] , 

m  (5, A)  =  sup  | q(x, A)  -  q(y, A) |  satisfies:  lim  m  (5,A)  =  0  uniformly  for 
q  |x-yi<5  SO  q 

A'  S  A  <  A".  Define  m  (5)  =  sup  m  (5, A).  Then  lim  m  (5)  =  0.  Since 

q  A'  <A<A"  q  SO  q 

-P-^  ’  ^  -  is  continuous,  there  is  a  constant  A  =  A(A',A"),  so  that 

m  (5, A)  <  AS,  for  A'  <  A  <  A".  Thus  m  (5)  <  AS. 

P  P 

Lemma  A. 2.  Let  A^  <  A'  <  A"  <  A^.  For  each  integer  k  >  0,  there  is  a 

function  w(5)  =  w^ ( 5 ) ,  defined  for  S  >  0,  such  that 

(1)  limw(5)  =0,  and 
SO 

(2)  |pR  h(A)  -  pk  h(A) |  <  w(h),  for  A'  <  A  <  A". 

Proof.  We  shall  use  the  Rayleigh  quotients 


n,,  ,  B  ( A ;  u , u ) 

R(  A,  u)  =  R(  A;  u)  =  - - — , 

<u,u>  <u,u>, 

h 


for  A'  S  A  <  A", 


u  e  S,  . 


The  eigenvalues  ^(A)  anc*  are  ^e^erm^ne<^  by 


(A. 10) 

Mk  hU) 

=  min  max  R( A; u) , 

A,  11 

dimU=k  ueU 

(A. 11) 

=  min  max  R(A;u), 

In.  >  11 

dimU=k  ueU 

where  U 

is  a  subspace  of  S^. 

We  may  assume  that 

(A. 12) 

<u,u>  =  1,  Hull2 
H 

<  M  and  |R(A;u) |  <  M 

(0,1] 

in  (A. 10).  Under  these  assumptions,  we  shall  estimate  |R(A; u)  -  R(A; u) | .  We 
shall  use  the  notations 


D  =  -  p(0,A)u(0)2  +  p(l,A)u(l)2, 


(A.  13) 


E  =  f  p(x, A)u' (x)2dx,  F  =  ["  q(x,  A)u(x)2dx, 

I  J  n  J  n 


|G  =  J  u(x)  dx. 

(Under  our  assumptions,  G  =  1,  but  we  shall  not  use  this  until  later.)  Let 
E  denote  the  midpoint  quadrature  of  E,  and  let  F,G  denote  the  trapezoid 
quadratures  of  F,G,  respectively.  Furthermore,  let 

(A. 14)  e  =  E-E,  f  =  F-F,  g  =  G  -  G. 

Note  that 


(A. 15) 


Thus  R( A; u)  = 


R( A; u)  = 


D+E-F 


R(  A;  u)  = 


D+E-F 


D+E-F  D+(E+e)-(F+f ) 


Using  the  Taylor  expansion, 


J_  =  I  n 

G+g  G 


G(  1+0) 


1 _ ] 

,„,2  ’ 


.V.V.V.V.V' "J 


where  0  lies  between  0  and  we  obtain 

R(A;u)  =  hi - — — 5]  [D+(E+e)-(F+f )  ] 

G  G( 1+0)2 


_  D+E-F  g 

G  g( 1+e)2 


[— t]  +  h  1  -  — 1 - p]  (e-f ) 

l-  G  J  G  C(l+0)2 


Using  (A. 12),  this  implies 


(A. 16) 


|R(A;u) -R(A;u)|  <  — +  1 1 - 

(i+e)  ( l+e) 


le  -  f  | 


where  0  lies  between  0  and  g.  We  shall  now  estimate  |e|,|f|,  and  |g|. 

A  function  u  e  S.  is  continuous  and  piecewise  linear  on  [0,1).  On  the 
h 

interval  [x^,x^  +  ^],  u(x)  =  u^  +  u'fx-x^,  where  =  uCx^)  and 

u'  =  (u^j-Ujl/h.  By  definition 

1  n-1  n-1 

E  =  J  p(x, A)u' (x)2dx  =  ^ (u' )2JXl  +  1p(x,  A)dx,  and  E  =  ^ (u')2p^  i(A)h. 

°  i=0  Xi  i=0  +2 

n-1 

Therefore  e  =  E-E  =  ^  (u' J2^^1  (p^  t(A)  -  p(x, A) )dx,  and 

1=0  Xi  2 

n-1  n-1  1 

I  e  |  <  m  (*)hY  (u()2.  Furthermore,  h  V'  (uj)2  =  f  u'(x)2dx  S  ||u||2  <  M 

p2  L*  1  ^  i  Jn  H1  [0,1] 


(by  (A. 12)),  and  nM^)  -  Therefore 


(A.  17) 


i  \  <  ma 
I  e  1  S  — h. 


1  n-1 

r  2  V-'  f^l  +  l  2 

We  shall  consider  |g|  next.  G  =  J  u  dx  =  [ u^+u' (x-x.  )  ]  dx  = 


i=0  Xi 


h  2  2  .  h  Vr  2  2>h3V^/\2  r  h  2^  2  \ 

3  Z.(ui+uiui+i+ui+i)  =  2  Z-'Wi1  "r  Z/V  •  G  =  2  Wl5 

i=0  1=0  i=0  1=0 

3  n_1 

Therefore  g  =  G-G  =  —  ^(u')2  -  |rh^.  Thus  we  have  shown 


Theorem  A. 1 .  Suppose  that  (1.1)  satisfies  the  monotonicity  assumptions 


(Ml)  -  (M4).  Then  for  each  eigenvalue  of  (1.1),  there  is  an  >  0, 

such  that  for  0  <  c  <  cQ,  and  for  small  h,  the  discrete  problem  (A. 4)  has 
a  unique  eigenvalue  A^  ^  in  the  interval  (A^-e,A^+e). 

Proof.  Let  <  A'  <  \<  A"  <  A2‘  Lemmas  A- 1  and  A' 2  imP1y  that 
(A. 22)  |jlk  h(A)  -Mk(A)l  ^  2k(h),  for  A'  <  A  £  A", 

where  zR(h)  =  wR(h)  + C(k)h2.  Let  z(h)  =  max(zk_1 (h) ,  zR(h) , zR+J (h) ) .  Then 
for  j  =  k-l,k,k+l, 

( A.  23 )  |pj  h(A)  -Pj(A)|  <  z(h),  for  A'  i  A  <  A". 

Note  that  lim  z(h)  =  0. 
h-*0 

Now  since  Ak  is  an  eigenvalue  of  (1.1),  pk(Ak)  =  Since  pk(A)  is 
a  simple  eigenvalue,  Pfc-l^k^  K  ^  <  ^k+l^k^'  ^et 

m  =  mln( |pkl(Ak) | ,pk+1(Ak) ),  and  choose  e0  >  0  small  enough  so  that 
I Pk_i (A) I  >  m/2  and  pk+1(A)  2:  m/2  for  -  c0  <  A  <  Ak + cQ.  Let 
0  <  c  <  Eq.  The  monotonicity  assumptions  (Ml)  -  (M4)  imply  that  Mk(A)  is 
strictly  decreasing.  Therefore  (ik(Ak~c)  >  0  >  pk(Ak+c).  Now  choose  h^ 
small  enough  so  that 

(A. 24)  z(h)  <  m/4,  for  0  <  h  <  hQ, 

(A. 25)  Pk(Ak-c)  -  z(h)  >0,  for  0  <  h  <  hQ, 

and 

(A. 26)  Pk(Ak+c)  + z(h)  <  0,  for  0  <  h  <  hQ. 

The  inequalities  (A. 23)  -  (A. 26)  imply  (for  0  <  h  <  h^) 

(A.27)  ;k.h<V‘>  >  i‘k.h'Vc)  <  “• 


(A.  28)  nk_i  ^(A)  <  T°r  ^k_G  5  S  A^+e, 

and 

(A. 29)  p,  ,  (A)  >  0,  for  A ,-e  <,  A  <  A ,+e. 

k+l,h  k  k 

This  implies  that  p^  ^(A)  has  a  zero  ^k  h  *n  (A^-e, A^+e) ,  anc*  T°r 
i  *  k,  p^  ^(A)  ^as  n°  zero  *n  ( A^-e, A^+e ) .  Q-E.D. 

Remark .  If  (1.1)  satisfies  the  monotonicity  assumptions  (Ml)  -  (M4),  then  the 
pure  finite  element  eigenvalue  p^  ^(A)  is  a  strictly  decreasing  function. 

Let  A  be  the  unique  zero  of  p  (A).  The  minimax  principle  implies  that 

iv  i  n  a  }  ii 

p^(A)  <  p^  ^(A),  and  therefore  A^  <  A^  These  facts  may  not  be  true  for 
Pk  ^(A)  and  A^  h-  IT  ^  \  h  is  always  true,  then  Theorem  4.3  is  also 

valid  in  the  case  that  A  is  an  eigenvalue  of  (1.1). 


Theorem  A. 2.  Let  A^  <  AQ  <  A^,  and  suppose  that  AQ  is  not  an  eigenvalue 
of  (1.1).  Then  there  is  an  e  >  0,  such  that  for  small  h,  the  discrete 
problem  has  no  eigenvalue  in  the  interval  ( AQ-e, Ag+e) . 


Proof.  If  the  statement  is  false,  then  there  is  a  sequence  h^ — >0  and 
eigenvalues  A^  ^  of  the  discrete  problem,  so  that  A^  ^  — >A^.  Let 


n  n 


n  n 


u  (x)  be  an  eigenfunction  corresponding  to  A,  (with  ||u  ||  . 

n  K.n  n  ..i  r  n  m 

n  n  H  10,  1 


=  1). 


Thus  B.  (A,  ,  ;u  ,v)  =  0,  for  all  v  e  S,  .  Since  the  u  are  bounded  in 

h  k  ,  h  n  h  n 

n  n  n  n 

H ^ [ 0 , 1 ) ,  there  is  a  subsequence  (which  we  again  denote  u  )  which  has  a  weak 


n 


.1. 


limit  u„  in  H  [0,1].  Since  h  — » 0  and  A.  — >An,  it  follows  that 

0  n  k  ,  h  0 

n  n 

B(A_;u„,v)  =  0  for  v  6  U  S,  .  Since  US,  is  dense  in  H*[0,1],  this 
(JO  h  n 

n  n 

implies  that  B(Aq;Uq,v)  =  0  for  all  v  e  H  ^ ( 0 , 1 ] .  But  this  means  that  A 


is  an  eigenvalue  of  (1.1). 


*  1 
► 

K 5L 


:i 


Q.E.  D. 


a 


v  w-  y- 
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